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ABSTRACT: 

This  paper  generalizes  the  classical  momentum  theory  as  usually  applied 
to  propellers, windmills  and  lifting  rotors  into  a  single  unified  treatment. 
It  also  extends  the  analysis  to  include  the  regime  in  which  flow  through 
part  or  all  of  the  actuator  disc  is  reversed  with  respect  to  the  remote  flow 
field.   Dimensional  analysis  is  used  in  a  systematic  manner  to  reduce  the 
final  results  to  their  simplest  and  most  significant  forms.   It  is  shown 
that  the  performance  of  these  devices  can  always  be  represented  by  a  single 
parameter  family  of  curves  in  which  the  parameter  expresses  the  extent  to 
which  the  performance  of  the  actual  device  falls  short  of  the  theoretical 
limit.   Detailed  algebraic  solutions  are  derived  in  closed  form  for  the 
performance  of  both  the  idealized  zero  torque  and  finite  torque  actuators; 
these  represent  performance  limits  which  any  comparable  real  device  may 
approach  but  never  exceed.  A  qualitative  analysis  is  presented  concerning 
the  development  of  the  vortex  ring  state.   The  fundamental  dynamic  stability 
of  the  ideal  lifting  rotor  is  also  analyzed. 

While  this  paper  deals  with  a  classical  topic  in  fundamental  fluid 
mechanics,  the  approach  employed  is  original  and  many  of  the  results  derived 
here  are  in  a  form  that  is  essentially  new.   In  particular,  the  solutions 
obtained  for  the  reversed  and  transition  flow  regimes  and  for  the  idealized 
rotating  actuator  are  believed  to  represent  new  formulations  of  these  basic 
phenomena. 
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1.   Objectives  of  the  Analysis 

The  simple  actuator  disc  has  long  served  as  a  useful  idealization  and 
standard  of  comparison  for  the  performance,  within  a  certain  restricted 
operating  range,  of  a  "broad  class  of  fluid -mechanical  devices.   These  include 
the  ordinary  aircraft  propeller,  the  lifting  rotor  such  as  the  helicopter 
rotor  under  certain  conditions  of  vertical  ascent  or  descent,  and  the 
unshrouded  axial  flow  turbine  or  "windmill"  as  it  is  termed  in  this  paper. 
The  fundamental  analysis  presented  here  is  applicable  to  all  three  of  these 
devices  which  we  designate  here  collectively  as  "actuators." 

The  chief  purpose  of  the  present  analysis  is  to  extend  the  elementary  and 
classical  momentum  theory  of  the  idealized  actuator  so  as  the  obtain  theoretical 
performance  limits  over  the  entire  possible  operating  range,  including  the 
regime  in  which  the  flow  through  part  or  all  of  the  actuator  disc  is  reversed 
with  respect  to  the  remote  flow  field!   In  all  cases  considered  in  this  paper, 
however,  the  axis  of  the  actuator  remains  parallel  to  the  undisturbed  flow  far 
away. 

Another  purpose  of  this  paper  is  to  utilize  dimensional  analysis  in  a 
systematic  manner  to  reduce  the  final  results  to  their  simplest  and  most 
significant  forms.   It  will  be  shown,  in  particular,  that  the  performance  of 
any  of  the  foregoing  devices  can  always  be  represented  by  a  single  parameter 
family  of  curves;  this  parameter  expresses  the  extent  to  which  performance  of 
the  actual  device  falls  short  of  the  limit  theoretically  attainable  from  the 
comparable  ideal  actuator . 

This  paper  also  considers  certain  fundamental  aspects  of  the  dynamic  and 
hydrodynamic  stability  of  the  ideal  actuator. 


2.  Summary  of  Principal  Concepts,  Equations  and  Generalized  Performance 

Charts 

The  main  results  of  this  analysis  are  summarized  in  subsections  2.1  through 
2.11  inclusive.   The  detailed  derivations  are  given  elsewhere  in  this  paper, 
but  enough  information  is  included  in  these  subsections  to  define  all  results 
without  the  need  to  refer  to  the  actual  derivations . 

2.1  Nomenclature  for  Section  2 

All  symbols  which  occur  in  the  summary  of  principal  results,  subsections 

2.2  through  2.11  inclusive,  are  defined  in  the  text  when  first  mentioned  and 
also,  for  ease  of  reference,  in  the  summary  below.   Symbols  are  listed  below 
according  to  the  subsection  and  in  the  approximate  order  in  which  they  are 
introduced  in  the  text.  A  dash  --  after  the  definition  signifies  that  the 
quantity  is  dimensionless . 

Subsection  2.2 

F  =  force  exerted  by  actuator  upon  fluid,  positive  if  in  the  same  sense 

as  the  velocity  of  the  undisturbed  fluid,   lbf 
L  =  reaction  exerted  by  the  fluid  upon  a  lifting  rotor,  always  upward  and 

positive  by  definition,   lbf 
V  =  axial  velocity  of  undisturbed  fluid  relative  to  actuator  disc ,  always 

non-negative  by  definition,   ft /sec 
H  =  rate  of  vertical  climb  of  a  lifting  rotor,  positive  upward,   ft /sec 
t  =  sign  changing  operator  defined  by  Eq.  (2.2.3).  -- 

Subsection  2.3 

o 
p  =  fluid  density  (constant),   slugs/ft 

D  =  diameter  of  actuator,   ft 

P.  =  shaft  power  of  an  ideal  actuator,  positive  if  supplied  to  actuator. 

ft  lbf /sec 

P  =  actual  shaft  power,  positive  if  supplied  to  actuator,  ft  lbf /sec 


P  =  power  loss  due  to  all  non-ideal  effects,  always  non-negative. 

ft  lbf/sec 

N  =  shaft  rotational  speed,  positive  by  definition,   rev/sec 

\i     =  fluid  viscosity,   lbf  sec /ft 

Subsection  2.U 

Note:  All  dimensionless  coefficients  listed  immediately  below  are  defined 

mathematically  in  Table  2.U.I. 

C_  =  conventional  force  coefficient  in  p,  N,  D  system. 

C   =  conventional  power  coefficient  in  p,  N,  D  system. 

J   =  conventional  advance  ratio  in  p,  N,  D  system. 

Cp  =  power  loss  coefficient  in  p,  N,  D  system, 
f 

C.-,  C-p,  C   ,  C    =  performance  coefficients  in  the  p,  D,  V  system. 
Cp_,  C_v,  C   ,  C„   =  performance  coefficients  in  the  p,  D,  |p|  system. 
C_p,  C  • ,  C   ,  C    =  performance  coefficients  in  the  p,  D,  L  system. 

Subsection  2 . 3 

V,  =  axial  velocity  at  actuator  disc,  positive  if  in  the  same  sense  as  the 
d 

undisturbed  flow,  ft /sec 
cp,  t,  o\  ,  ct  =  sign  changing  operators  defined  by  Eqs.  (2.5). 
Subsection  2.6 
P,  W,  T,  R  =  symbols  denoting  propulsive,  windmill,  transition  and 

reversed  flow  operating  modes,  respectively. 
M   =  symbol  denoting  operating  point  which  corresponds  to  generation  of 

maximum  possible  windmill  power  |C,p|. 
B   =  symbol  denoting  operating  point  which  corresponds  to  maximum  drag 

force  at  zero  net  shaft  power. 
a,m  =  parameters  occurring  in  the  cubic  solution  for  thrust  coefficient 

as  defined  in  Eq.  (2.6.U)  of  Table  2.6.2.   -- 


Subsection  2.7 

TR,  TW  =  symbols  denoting  subregions  of  the  transition  mode  corresponding, 

respectively,  to  positive  and  negative  net  shaft  power. 

a,   m  =  parameters  occurring  in  the  cubic  solution  for  thrust  coefficient 

as  defined  in  Eq.  (2.7-7)  of  Table  2.7-H-  -.-. 

B  =  parameter  occurring  in  the  cubic  solution  for  velocity  coefficient 

as  defined  in  Eq.  (2.7-7)  of  Table  2.7-i+-   -- 

Subsection  2.8 

S  =  symbol  denoting  operating  point  which  corresponds  to  static  conditions 

(V=  0). 

8  -   parameter  occurring  in  the  cubic  solution  for  rate  of  descent  at  zero 

power  as  defined  in  Note  2  of  Table  2.8.1. 

Subsection  2.9 

r  =  radius  of  non-swirling  core  in  remote  slipstream,   ft 

cj\    ,  C_\    =  maximum  values  of  conventional  thrust  and  power  coefficients 
J  MAX  Jmx 

that  can  be  attained  at  any  specified  advance  ratio  J  by 

an  idealized  finite  torque  actuator. 

R  =  slipstream  radius,   ft 

r  =  slipstream  radial  coordinate,   ft 

U)  =  shaft  angular  velocity,   rad/sec 

V  =  slipstream  axial  velocity,   ft /sec 

V  =  slipstream  peripheral  velocity,   ft /sec 

T  =  slipstream  vortex  strength  outside  core  region,   ft  /sec 

k=  slipstream  advance  ratio  defined  by  Eq.  (2.9-2)  of  Table  2.9-1. 
A  =  slipstream  loading  parameter  defined  by  Eq.  (2.9-3)  of  Table  2.9-1- 

V  =  slipstream  swirl  parameter  defined  by  Eq.  (2.9.U)  of  Table  2. 9.1. 
TL  =  relative  radius  of  non-swirling  core  as  defined  by  Eq.  (2.9- 5)  of 

Table  2. 9.1.   — 
C  =  auxiliary  function  defined  by  Eq.  (2.9.7)  of  Table  2. 9.1.   -- 

k 


Subsection  2.10 

6L  =  perturbation  in  lift,   lbf 
W  =  total  weight  of  rotorcraft.   lbf 
t  =  time  coordinate,   sec 
9  =  dimensionless  time  coordinate  defined  by  Eq.  (2.10.U)  of  Table 

2.10.1.   -- 
C   ,  C  •  ,   C_  =  dimensionless  performance  coefficients  in  p,  D,  W 

system  as  defined  by  Eqs.  (2.10.1)  through  (2.10.3)  of 
Table  2.10.1.   — 
X,  Y  =  dimensionless  stability  derivatives  defined  by  Eqs  (2.10.5)  through 

(2.10.8)  of  Table  2.10.1.   — 
(j,,  v  =  perturbation  acceleration  and  velocity  response  constants,  respec- 
tively, of  control  system  as  defined  by  Eq.  (2.10.9)  of  Table  2.10.1. 
3  =  exponential  response  rate  of  lifting  rotor  to  a  small  perturbation. 
(Operation  is  stable  if  B  is  negative). 
2.2  The  Lifting  Rotor  in  Vertical  Ascent  or  Descent 

In  dealing  with  the  lifting  rotor  in  vertical  ascent  or  descent,  the 
rotor  axis  is  placed  vertically,  and  the  sense  of  the  flow  through  the  actuator 
is  so  chosen  that  the  reaction  force  F  exerted  by  the  fluid  upon  the  actuator 
itself,  whether  it  be  a  propulsive  thrust, a  windmilling  drag  or  a  power-on 
reversed  thrust,  always  acts  upward  and  thereby  constitutes  a  useful  lift  L. 
These  circumstances  in  no  way  affect  the  flow  pattern  itself,  nor  the  basic 
method  of  analysis. 

It  is  convenient,  however,  to  introduce  two  minor  changes  of  variable 
when  dealing  specifically  with  the  lifting  rotor.   In  place  of  the  usual  axial 
force  F  which  may  be  either  positive  or  negative,  we  introduce  the  lift  L  which 
is  always  directed  upward,  and  which  we  therefore  take  as  positive  by  definition. 


In  place  of  the  remote  velocity  V  which  is  always  non-negative  by  definition, 
we  introduce  the  rate  of  climb  H  which  is  positive  for  vertical  ascent,  negative 
for  vertical  descent.   (The  case  of  slanting  ascent  or  descent  falls  outside 
the  scope  of  this  analysis.) 

The  required  changes  of  variable  are  simple  accomplished  by  the  substitutions 

F  =  t  L  (2.2.1) 

V  =  t  H  (2.2.2) 

where  t  is  a  sign  changing  operator  defined  as  follows,  namely, 

T  =  -[fy  =  ±i  (2.2.3) 

By  means  of  these  substitutions ,  any  result  originally  derived  for  the 
propeller  or  windmill  can  be  converted  to  the  corresponding  result  for  the 
lifting  rotor  in  vertical  ascent  or  descent. 
2.3  Real  versus  Ideal  Actuators 

The  present  analysis  is  based  on  comparing  the  performance  of  any  real 
propeller,  windmill  or  lifting  rotor  with  that  of  a  corresponding  idealized 
zero  torque  actuator.  As  illustrated  in  Section  11  of  this  paper,  the  above 
comparison  can  be  made  in  more  than  one  way.  However,  the  scheme  outlined 
below  appears  to  be  by  far  the  simplest  and  most  versatile.   Hence  this  is 
the  basic  method  which  is  consistently  applied  throughout  this  paper,  except 
for  the  auxiliary  discussion  in  section  11. 

In  this  method,  the  idealized  actuator  is  arbitrarily  assigned  values  of 
p,  D,  V  and  F  which  are  respectively  equal  to  those  of  the  real  device,  by 
definition.   This  then  suffices  to  fix  the  theoretical  shaft  power  P.  of 
the  ideal  actuator  according  to  a  definite  relation  which  is  derived  elsewhere 
in  this  paper  and  which  can  here  be  symbolized  in  the  form 


P.  =  P.  (p,  D,  V,  F)  (2.3.1) 

On  this  basis,  the  shaft  power  P  of  the  real  actuator,  whether  positive 
or  negative,  is  always  algebraically  greater  than  the  ideal  shaft  power  P.. 

The  difference  between  these  two  quantities,  which  we  denote  by  symbol 
P  ,  represents  a  relative  loss  of  power  due  to  all  non-ideal  effects,  including 
those  associated  with  slipstream  rotation,  viscosity,  number  of  blades,  detailed 
blade  geometry,  flow  unsteadiness,  separation  effects  and  so  on.  Hence  P„  is 
in  general  a  complex  function  not  only  of  p,  D,  V  and  F,  but  also  of  shaft 
speed  N,  of  viscosity  \i,   and  of  various  geometric  parameters.   It  is  necessarily 
non-negative.   Thus  symbolically, 

Pf  =  Pf  (p,  V,  D,  F,  N,  n, )  *  0  (2.3.2) 

Therefore  the  shaft  power  P  of  the  real  actuator  may  finally  be  represented 
in  the  form 


P  =  P±  +  Pf  (2.3.3) 


Since  the  principal  aim  of  this  paper  is  to  establish  a  simple  performance 
limit  which  any  real  actuator  can  approach  but  never  exceed,  attention  is^ 
focussed  primarily  on  the  elementary  analytical  relation  which  is  symbolized 
by  Eq.  (2.3.I)  and  which  represents  the  idealized  zero  torque  actuator  itself. 
The  more  complex  and  obscure  relation  symbolized  by  Eq.  (2.3-2),  which  relates 
the  real  actuator  to  the  ideal,  is  defined  conceptually  in  this  paper,  but  is 
only  partially  analyzed.   In  a  later  section  the  overall  power  loss  P„  is 
broken  down  further  into  a  rotation  loss  P„  plus  a  residual  loss  P  .   Of  course, 
the  latter  quantity  is  in  principle  amenable  to  further  analysis,  but  this 
requires  resort  to  methods  which  are  either  semi-empirical  in  character,  or 


else  far  more  complex  than  those  considered  in  this  paper. 

In  general,  for  any  well  designed  actuator  operating  at  or  near  its 
design  point,  the  overall  power  loss  P  should  be  comparatively  small  and  well 
suited  to  estimation  in  any  specific  instance  by  accepted  engineering  methods 
based  on  the  utilization  of  appropriate  test  data. 

Attention  is  also  invited  to  the  role  of  shaft  speed  N  in  the  situation. 
Note  that  power  loss  Pf  is  a  function  of  N,  as  indicated  in  Eq.  (2.3.2),  but 
that  the  basic  relation  (2.3.1)  for  the  ideal  power  P.  itself  is  actually 
independent  of  N!   This  fact  is  subsequently  used  to  simplify  the  basic  analysis 
in  a  very  significant  way. 
2.k     Dimensionless  Performance  Coefficients 

The  analysis  of  any  propeller,  windmill  or  lifting  rotor  can  be  simplified 
and  generalized  by  reducing  all  quantities  and  equations  to  non-dimensional  form. 
Since  the  three  fundamental  dimensions  of  force,  length  and  time  are  involved 
in  this  problem,  it  is  necessary  for  this  purpose  to  choose  three  of  the  basic 
parameters  of  the  problem  as  dimensional  reference  parameters.   This  can  be  done 
in  various  ways.   Four  of  the  most  important  possibilities  are  summarized  in 
Table  2.1+.1. 

By  far  the  most  commonly  used  system  of  dimensionless  coefficients  is  that 
labelled  the  "Conventional  System"  in  the  table.   This  uses  p,  N  and  D  as 
reference  parameters  and  results  in  the  three  familiar  coefficients  commonly 
known  as  the  thrust  coefficient  C„,  power  coefficient  C  and  advance  ratio  J, 
as  defined  in  the  table. 

We  observe,  however,  that  while  shaft  speed  N  is  a  significant  parameter 
of  the  overall  problem,  it  is  not  involved  at  all  in  the  fundamental  relation 
of  the  idealized  zero  torque  actuator  as  symbolized  by  Eq. (2.3.1).   Hence, 
although  this  relation  can  certainly  be  non-dimensionalized  in  the  conventional 
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p,  N,  D  system,  the  results  so  obtained  will  not  be  expressed  in  terms  of  the 
fewest  possible  number  of  dimensionless  parameters.   In  effect,  the  conventional 
system  introduces  an  extraneous  parameter  N  into  the  formulation  of  the  basic 
relation  (2.3.1)  and  thereby  creates  one  superfluous  degree  of  freedom. 

The  above  complication  is  undesirable  and  can  be  avoided  simply  by  dropping 
the  conventional  p,  N,  D  system  and  adopting  instead  any  one  of  the  three  alter- 
native schemes  shown  in  Table  2.U.I.   Notice  that  the  various  dimensionless 
coefficients  that  occur  in  these  alternative  systems  can  all  be  expressed  very 
simply  in  terms  of  the  familiar  coefficients  CL,,  C  and  J  of  the  conventional 
system.   These  simple  conversions  are  summarized  in  the  table. 

The  above  three  basic  systems  are  theoretically  equivalent  in  that  they 
all  express  the  same  underlying  relationship  in  terms  of  the  same  number  of 
dimensionless  parameters.  The  respective  detailed  results,  however,  are  not 
necessarily  equally  apt  and  convenient  in  every  possible  application.  As  a 
rule,  the  specific  context  determines  which  one  of  these  three  systems  expresses 
the  pertinent  results  in  the  most  appropriate  form. 

It  is  also  worth  noting  that  under  any  circumstances  in  which  parameter  N 
reenters  the  problem  in  a  significant  way,  the  conventional  coefficients  CL,, 
Cp  and  J  once  again  become  relevant.  This  is  true,  for  example,  of  the  idealized 
finite  torque  actuator  as  analyzed  in  detail  in  section  9,   and  summarized  in 
subsection  2.9. 
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2.5  Sign  Changing  Operators 

Owing  to  the  extremely  wide  range  of  operating  conditions  considered  in 
this  paper,  the  velocity  V,  through  the  actuator  disc,  the  axial  force  F,  the 
ideal  shaft  power  P.  and  actual  shaft  power  P  are  each  subject  to  changes  in 
algebraic  sign.  As  an  aid  in  dealing  with  these  changes  in  an  orderly  fashion, 
we  introduce  a  group  of  sign  changing  operators  defined  as  follows: 


T--|5p+I  (2.5.2) 


P. 

-nrr=  +  i  (2.5.3) 


i 


K 


"-•J5J--  +  1  (2.5.10 

The  usefulness  of  these  operators  for  generalizing  the  basic  analytical 
relations  will  become  apparent  in  the  following  subsections . 
2.6  The  Basic  Solution  in  the  p,  D.  V  System 

The  basic  momentum  and  energy  analysis  shows  that  the  complete  operating 
envelope  of  an  idealized  zero  torque  actuator  may  be  subdivided  into  four 
distinct  ranges  as  indicated  in  Table  2.6.1.  Each  operating  range  has  its  own 
characteristic  idealized  flow  pattern  as  shown  schematically  in  Fig.  2.6.1. 
The  two  basic  equations  presented  in  Table  2.6.1  suffice  to  define  the  performance 
over  this  entire  range.   Provided  only  that  the  proper  respective  sign  changing 
operators  cp  and  t.  be  used  as  defined  in  the  table,  the  single  generalized 
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Table  2.6.1  Basic  Actuator  Performance  Relations 

As  Functions  of  Force  Coefficient  in  p,  D,  V  System 


Operating 
Mode 

Symbol 

Range 

Sign  Changing 
Operators 

Basic 
Eq 

<P 

a. 

i 

Propuls  ive 

P 

o  ^c1F  *  +  =0 

+1 

+1 

(2.6.1) 

Windmill 

W 

-  g  o/l  -  D  *cir,*o 

+1 

-1 

(2.6.1) 

Transition 

T 

-g(V3  +  D  ^c^s-JW-i) 

Not 
Applicable 

(2.6.2) 

Reverse 

R 

"   "^   S-gG/I+D 

-1 

+1 

(2.6.1) 

Basic  Equations 

Modes   P,  W  and  R 

2.6.1) 

V  ■(¥)  ♦•!  At)  +4  ciF3  +  ^  (: 

Mode  T 

f       -          n      _  p        ip                   fo^oN 

lp  "       3/3"        1F        1Pf 

Note:  For  definitions  of  coefficients,  see  Table  2.U.I. 
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Eq.  (2.5.1)  applies  to  all  three  ranges  designated  as  modes  P,  W  and  R!  The 
separate  relation  shown  in  Eq.  (2.6.2)  then  applies  over  the  remaining 
transition  range,  mode  T. 

Eq.  (2.6.1)  is  completely  consistent  with  the  usual  elementary  momentum 
analysis  of  the  simple  propeller  or  windmill.   Hence  it  is  not  surprising  that 
the  same  basic  result  is  applicable  to  both.  But,  as  Fig.  2.6.1  plainly  shows, 
the  range  of  fully  reversed  flow,  mode  R,  involves  a  flow  pattern  which  seems 
very  much  more  complex  than  the  simple  patterns  associated  with  modes  P  or  W. 
It  is  therefore  perhaps  surprising  that  Eq.  (2.6.1)  can  include  all  three  of 
these  modes  within  a  single  generalized  expression.   In  addition  to  generalizing 
the  classical  theory  in  this  striking  way,  the  present  analysis  also  yields  the 
novel  result  for  the  transition  range  as  represented  by  Eq.  (2.6.2);  this 
particular  relation  has  no  counterpart  in  the  usual  elementary  theory! 

The  basic  overall  relation  defined  jointly  by  Eqs .  (2.6.1)  and  (2.6.2)  is 

shown  in  Fig.  2.6.2  as  a  family  of  curves  of  C.  _.  versus  C   ,  with  C    as  para- 

lir  _Lb         -LrT 

meter.  The  curve  for  C    =  0  represents  the  ideal  actuator.  This  defines  the 
limiting  performance  which  any  real  actuator  can  approach  but  never  exceed. 

The  loss  parameter  C.  __  expresses  quantitatively  how  much  the  performance 
of  any  real  device  falls  short  of  the  theoretical  ideal.   It  also  provides  a 
convenient  basis  for  comparing  the  relative  performances  of  various  real  actuators 

Two  critical  points  on  the  ideal  curve  are  worthy  of  note.   Point  M  shows 

that  the  magnitude  of  the  power  coefficient  |cip|  generated  by  any  windmill  can 

2n 
never  exceed  g=r  =  0.233.   Point  B  shows  that  the  magnitude  of  the  drag  coefficient 

|C, „  \   which  can  be  generated  without  the  expenditure  of  net  shaft  power  can 

never  exceed  — 7=  =  O.605.   For  a  real  actuator  of  specified  constant  C-ipf>  these 
3v  3 

magnitudes  decrease  to  (~  -  C-iPf.)  and  ( — —  -  C.^),  respectively. 
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It  should  "be  mentioned  that  if  we  set  C1pf,  =  0  in  Eq.  (2.6.1),  this 
analytical  expression  can  generate  not  only  the  major  part  of  the  ideal  per- 
formance envelope  shown  in  Fig.  2.6.2,  but  also  certain  extraneous  solutions 
which  if  they  were  shown  here  would  fall  within  the  interior  region  of  the  chart. 
These  details  are  in  face  shown  later  in  Fig.  6.2.   It  is  demonstrated  elsewhere 
in  this  paper  that  one  of  these  branches,  although  it  formally  satisfies  the 
equation,  represents  an  unstable  solution  which  cannot  be  physically  realized 
in  the  present  context.   The  other  extraneous  branches  are  theoretically  stable 
but  represent  performance  which  is  less  than  optimal.   These  extraneous  solutions 
are  not  considered  further  in  this  summary  of  principal  results. 

There  is  one  singular  operating  condition,  however,  which  does  warrant  brief 

comment  at  this  point.   That  is  the  so-called  ring  vortex  state,  represented  in 

Fig.  2.6.2  by  point  V.   This  condition  is  of  particular  interest  in  connection 

with  the  lifting  rotor  because  it  is  associated  with  a  catastropic  loss  of  lift! 

The  typical  flow  pattern  of  the  vortex  ring  state  is  shown  schematically  in  Fig. 

2.6.3.  This  operating  state  is  considered  further  in  the  detailed  analysis. 

The  original  solution  summarized  in  Table  2.6.1  is  based  on  CnTn  as  the 
&  IF 

independent  variable;  this  choice  reduces  the  analytical  expressions  to  their 
simplest  form.   In  practice,  however,  it  is  frequently  required  to  take  CLp 
as  the  independent  variable.  Fortunately,  the  basic  equations  can  in  fact  be 
solved  for  C._  as  a  function  of  Cn_,  even  though  Eq.  (2.6.1)  happens  to  be  a 

.Lb  Xr 

cubic  equation  in  the  variable  C  .  When  extraneous  branches  of  the  solution 
are  discarded,  the  final  results  of  this  process  can  be  summarized  as  shown  in 
Table  2.6.2.  Note  that  the  single  equation  (2.6.1)  is  now  replaced  by  the  two 
expressions  (2.6.3)  and  (2.6.U)  which  are  applicable  over  different  ranges. 
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ACTUATOR  DISC 


APPROXIMATE   STREAMLINES 


Fig.  2.6.3    SCHEMATIC  DIAGRAM  OF  LIFTING  ROTOR 
IN  VORTEX  RING  STATE. 
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Table  2.6.2  Basic  Actuator  Performance  Relations 

as  Functions  of  Power  Coefficient  in  p,  D,  V  System 


Mode 


Symbol 


Range 


<P      Eq, 


Propulsive 


0  <.   (C 


IP 


ciPf >  *  + 


+l  ;  (2-6.3: 


Windmill 


W 


—  £  (c 

27   l  IP 


C    )  ^  0 
lPf 


N.A.  i  (2.6.U) 


Transition 


t(M*^-^*:zm\ «. 


(2.6.5) 


Reverse 


R 


tTO^-°,(t 


'IP    lPf' 


-1 


(2.6.3) 


Basic  Equations 


Modes  P  and  R 


C1F  ~  1?  (C1P  "  ClPf )  I 


1/3 


I 


27  1°1P   UlPf; 


-1 


1  -  A +  If  (cip  -  ciPf )"1 


1/3 
l/3l 


(2.6.3) 


Mode  ¥ 


=  /^¥ 


C1F  :      t  (C1P  "  ClPf  >  COS 


ry 

2tt1 

— 

+ 

m 

a 

3 

(2.6.U) 


i  where  cos  a  =  J-   p—  (C   -  C   )   ,  0  <  a   <:  ^  ,  and  where  m  has  the  values  -1 


and  +1,  respectively,  over  the  following  sub-ranges: 


(Table  continues  on  next  page) 
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Table  2.6.2          Cont'd 

m  = 

-1 

-  if  *  (cip  -  cipf }  *  ° 

m  = 

+1 

.5Wr^i(c     -  c     )  ^  -  2tt 
%[  j$  J       ip      iPf;        27 

Mode  T 

°1F  =  "  3v5  "  (°1P  "  ClPf ) 

(2.6.5) 

Note:  For  definitions  of  coefficients,  see  Table  2.U.I. 
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However,  the  analytical  solutions  summarized  in  Tables  2.6.1  and  2.6.2  are 

entirely  equivalent,  and  "both  generate  precisely  the  same  family  of  curves, 

namely,  that  family  previously  shown  in  Fig.  2.6.2. 

It  should  be  appreciated  that  the  form  of  solution  based  on  the  p,  D,  V 

dimensional  system  as  illustrated  in  Fig. 2.6.2  can  be  useful  in  any  situation 

in  which  V  is  both  known  and  non-vanishing.   The  windmill  is  perhaps  the  best 

example,  because  no  windmill  can  generate  useful  shaft  power,  nor  even  a  useful 

drag  force,  unless  V  is  non-zero.  Conversely,  this  particular  reference  system 

is  virtually  useless  for  dealing  with  any  situation  in  which  performance  under 

static  conditions,  V  =  0,  is  a  significant  factor  because  under  these  conditions 

all  four  of  the  dimensionless  coefficients  CnTP,  Cn_,  C    and  C   become  infinite 

.Lt1    J-Jr    _L.Pi       _L1M 

and  comparatively  meaningless.   The  important  cases  of  the  propeller  under  static 
thrust  or  the  lifting  rotor  in  stationary  hovering  flight  are  examples  of  situa- 
tions which  cannot  be  usefully  treated  within  the  framework  of  the  p,  D,  V 
dimensional  system.   To  deal  with  these  problems  we  must  employ  the  other 
alternative  systems . 
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2.7  The  Basic  Solution  in  the  p,  D,  1p|  System 

In  the  previous  p,  D,  V  dimensional  system,  the  parameters  p  and  D  are 
innately  positive,  and  V  is  restricted  to  non-negative  values  by  definition. 
We  note,  however,  that  in  the  present  dimensional  system,  shaft  power  P  itself 
may  be  either  positive  or  negative,  depending  on  the  particular  mode  of  operation 
considered.   It  is  essential  that  all  dimensional  reference  quantities  be 
restricted  to  positive  values  only,  hence  we  take  |pj  rather  than  P  as  the 
appropriate  reference  quantity. 

We  also  note  that  the  overall  transition  range,  mode  T,  must  now  be  sub- 
divided into  two  sub -intervals  corresponding,  respectively,  to  positive  and 
negative  values  of  net  shaft  power  P.  We  denote  these  as  modes  TR  and  TW, 
respectively. 

Upon  changing  from  the  p,  D,  V  system  to  the  p,  D,  |p|  system,  the  original 
pair  of  fundamental  equations  as  previously  given  in  Table  2.6.1  are  thereby 
converted  into  the  corresponding  pair  of  equations  now  shown  in  Table  2.7.1. 
Of  course  the  change  involved  is  purely  a  matter  of  algebraic  form;  the  under- 
lying physical  relationships  are  still  precisely  the  same  as  before. 

We  observe  that  Eq.  (2.7.1)  of  the  table  is  linear  in  Cp  but  cubic  in 

Cp„,  whereas  Eq. (2.7.2)  is  linear  in  Cp„  but  cubic  in  Cp  .  Each  equation  is 

here  solved  for  the  linear  term  thereby  expressing  the  final  relations  in  their 

simplest  algebraic  forms,  as  shown.  However,  it  is  sometimes  desirable  to  write 

both  these  relations  for  CAT1  in  the  form  of  an  explicit  function  of  C_T7..   In 

2F  2V 

this  case  it  becomes  useful  to  distinguish  between  the  power  input  and  power 
output  modes.  The  solutions  for  these  two  cases  are  summarized  separately  in 
Tables  2.7.2  and  2. 7-3- 

Conversely,  it  is  sometimes  convenient  to  be  able  to  solve  all  relations 
consistently  for  COT7.  as  a  function  of  C0„.  This  requires  the  inversion  of  Eq. 
(2.7.2)  which  is  summarized  in  Table  2.7.^. 
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Table  2.7.1  Basic  Actuator  Performance  Relations 


in  p,  D,  1p|  System 


Mode 


Symbol 


Range 


cp 


Power  Input  Modes  - 


Propulsive 


-  (1-C    Y 
2  U  U2Pf ; 


1/3 


*  C2F  *  ° 


+1 


+1 


(2.7, 


Reverse 


-  fl-C    V 
2  U  U2Pfj 


1/3 


'2F 


§  (^+i)(i-c2pfy 


1/3 


+1 


Transition 
(Power  In) 


TR 


I  (£)^r 


1/3 


£  C2y  £  +  - 


+1 


N.A. 


Power  Output  Modes 


Windmill 


W 


0  *  C2F  *  " 


I  (V5-D  d^2Pf)2 


1/3 


-1 


+1 


Transition 
(Power  Out) 


•  TW 


l(^)(1^ 


1/3 


*  C2V  *  +  " 


-1 


N.A, 


Basic  Equations  - 


Modes   P,  R,  W 


(1   -   oC2pf) 


'2V 


oC 


2F 


—  2F 

TTCP    (1     -     0C2pf) 


Modes   TR  and  TW 


llLi       2      a  (i  -  oC2pf) 
2F  "  "  3/3    2V    "  C2V 


Note:   For  definitions  of  coefficients,  see  Table  2 .U.l. 
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Table  2.7.2     Basic  Actuator  Performance  Relations   in  Terms   of 
Velocity  Coefficient  for  Power  Input  Modes   in   p,  D,    1p|  System 


Mode 


Symbol 


Range 


<P 


Eq. 


Propuls  ive 


0   *  C2V  *  +   " 


+1 


(2.7.3) 


Reverse 


R 


0    *  C2V   *  + 


££  (1   -  C         ) 

2tt  U         2Pf; 


1/3 


-1 


(2.7.3) 


Transition 
(Power  In) 


TR 


—  fl  -  C        ) 
2tt  U         2Pf ; 


1/3 


* C2V  *  +   " 


N.A. 


(2.7.10 


Basic  Equations 


Modes  P  and  R 


'2F  ~  ]£ 


?  (i  -  c2pf  )2r 


\l/3 


1  + 


1    +  2ttcp         °2V 

27     (1   -  C2pf  )j 


1/3 


•[- 


i    +  2ncp  2V 

27    (1  -  C2pf  Jj 


1/3 


(2.7.3) 


Mode  TR 


'2F  "  "  3V§"     2V 


a.     2         '  c2Pf ^ 


'2V 


(2.7.U) 


Note:     For  definitions  of  coefficients,   see  Table  2.U.1, 
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Table  2.7-3     Basic  Actuator   Performance  Relations   in  Terms   of 
Velocity  Coefficient  for  Power  Output  Modes   in    p,   D,     \~p\  System 


Mode 


Symbol 


Range 


Eq, 


Windmill 


W 


E  (1   +  C2Pf ) 


1/3 


*  C2V  *  +   * 


(2.7.5 


Transition 
(Power  Out) 


TW 


;f(^-^P^ 


£  + 


(2.7.6) 


Basic  Equations    - 


Mode  W 


'2F 


a 


f  (1   +  °2Pf '   °2V  COS 


'a  ,2tk 

3  +  m   (-) 


(2.7.5) 


/       (1  +  C        ) 
where   cos   a  =  Jt »s ?   °   ^  ^  ^  o"  j   and  where  m  has   the   values   -1 

U2V 

and  +1,  respectively,  over  the  following  sub-ranges: 


m  = 

-1 

'%  <1+C2pP 

1/3 
_    *  C2V  *  +  °° 

m  = 

+1 

[i  (1+c2Pf )" 

1/3   ^  C        * 

[;@£M 

1/3 

Mode  TW 


p   (l  +  C   ) 

c     =  _  -H-  c       +  2pf 

C2F     3^C2V        C2V 


(2.7.6) 


Note:  For  definitions  of  coefficients,  see  Table  2.U.I. 
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Table  2.7.^-  Basic  Actuator  Performance  Relation  in  Terms  of 
Force  Coefficient  for  Transition  Modes  in  p,  D,  |p|  System 


Mode 

Symbol 

Range 

Q1 

Eq. 

Transition 
(Power  In) 

TR 

- 

~£  (S  *  i)  (i  -  c2pf  f 

1/3  ■>-  - 

+1 

(2.7.7) 

Transition 
(Power  Out) 

TW 

-  h  ts.  - 1)  (i  +  c2Pf)2 

1/3 -2F*-- 

-1 

(2.7.7) 

Basic  Equation   - 

Modes  TR  and 
where 

TW 

CO 

■! 

3/2 

(2.7.7) 

C2V  *  2  ^  ("  C2F)   CO! 
f~  a  (1  -  o< 

s  >--|/-i . 

Note:   For  definition  of  coefficients,  see  Table  2.U.I. 
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The  form  of  the  solution  is  shown  graphically  in  Fig  2.7.1  for  the  power 
input  modes  and  in  Fig  2.7.2  for  the  power  output  modes.   In  each  case  the 
solution  consists  of  a  family  of  curves  of  C0  versus  C   with  the  power  loss 
coefficient  Cpp_  as  parameter.   Of  course  the  curves  for  Cpp_  =  0  once  again 
represent  a  limit  that  any  real  actuator  may  approach  hut  never  exceed. 

The  upper  part  of  Fig.  2.7.1  represents  a  propeller  in  normal  operation. 
The  lower  part  represents  a  propeller  under  conditions  of  power -on  reversed 
thrust.   Fig.  2.7-1  represents  a  windmill  generating  useful  power  and  undergoing 
a  drag  force.  We  can  think  of  this  case  as  depicting  a  windmill  of  fixed  size 
generating  fixed  shaft  power  under  various  wind  velocities.   Note  the  existence 
of  a  limiting  velocity  below  which  the  specified  power  cannot  he  attained.   This 


1/3 

.  Also  at  any  given 


[?7  P 

_—  (l  +  Cpp„) 
- 
value  of  Cpv  above  the  critical,  there  are  of  course  two  solutions  corresponding 

to  low  and  high  drag,  respectively. 
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VELOCITY  COEFFICIENT   C2p  = 
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Fig.    2.7.2    BASIC  ACTUATOR  PERFORMANCE  CURVES  FOR  POWER 
OUT  PUT  MODES  IN  />,  D,  |P|  SYSTEM. 
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2.8  The  Basic  Solution  for  the  Lifting  Rotor  in  the  p,  D,  L  System 

Straightforward  algebraic  conversion  of  the  original  pair  of  fundamental 
equations  of  Table  2.6.1  into  the  p,  D,  L  dimensional  system  produces  the 
corresponding  pair  of  equations  for  the  lifting  rotor  as  shown  in  Table  2.8.1. 
These  equations  show  power  coefficient  C~p  as  an  explicite  function  of  rate  of 
climb  coefficient  C  •  ,  with  power  loss  coefficient  CL   as  parameter.  This 
arrangement  not  only  yields  the  simplest  algebraic  result,  but  is  also  the 
most  convenient  for  various  other  purposes  as  well.   Of  course,  if  necessary, 
both  of  these  relations  can  be  readily  inverted  to  give  CL  •  as  an  explicit 
function  of  (C,.  -  C_  _);  however,  this  again  entails  solving  a  cubic 
equation  and  the  detailed  results  obtainable  in  this  way  are  not  included  here. 

The  fundamental  relation  defined  in  Table  2.8.1  is  shown  graphically  in 
Fig.  2.8.1  as  a  family  of  curves.   Once  again,  the  curve  for  C  ^  =  0  represents 
a  performance  limit  that  any  real  lifting  rotor  can  approach  but  never  exceed. 
Of  particular  interest  are  the  two  critical  points  marked  S  and  B.   Point  S 
shows  that  the  power  coefficient  C_p  required  to  maintain  steady  hovering  can 
never  be  less  than  J  —  =*  0.798-  Point  B  shows  that  the  rate  of  descent  at  zero  net 

TT  I 

shaft  power  as  measured  by  |C^J  can  never  be  smaller  than  J  -X— =  1.286.  The 
changes  of  these  values  with  increasing  loss  coefficient  C_  „  is  also  plainly  shown. 

It.  is  noteworthy  that  these  results  represent  an  elementary  solution  in 
closed  form  over  the  entire  operating  range!  Hitherto,  elementary  solutions 
based  on  the  ordinary  momentum  analysis  have  been  available  over  those  portions 
of  the  range  which  correspond,  in  the  present  terminology,  to  modes  P  and  W 
only.   Moreover,  the  present  results  embody  the  useful  parameter  CUpf  which 
facilitates  the  comparison  of  any  real  lifting  rotor  with  the  ideal,  or  with 
some  other  real  rotor. 
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Table  2.8.1  Basic  Lifting  Rotor  Performance  Relations  in  Terms 
of  Rate  of  Climb  Coefficient  in  p,  D,  L  System 


Mode 

Symbol 

Range 

<p 

T 

a. 

i 

Eq. 

Propulsive 

P 

+  »s  C  •  >  0 

5n 

+l 

+1 

+1 

(2.8.1; 

Reverse 

R 

oaC3H£-f£(^+l)}4 

-l 

-1 

+1 

(2.8.1 

Transition 

T 

-^(^i)j-^c3fi,-^(^-i)j-* 

N.A. 

N.A. 

N.A. 

(2.8.2 

Windmill 

W 

-(?(^-i)j-^c3]1,-- 

+1 

-1 

-1 

(2.8.1: 

Basic  Equations  - 

Modes  P,  R, 

W 

c 

(2.8.1J 

^--c-¥^JC-¥)i^y^ 

Mode  T 

TT       3 

3P    3^  3H   L3H   °3Pf 

(2.8.2 

Note  1:  At  C  •  =  0,  Eq.  (2.8.1)  gives  C   =  v|  +  c3Pf' 


Note  2:  At  C   =  0,  Eq.  (2.8.2)  gives  C  • 


M 


1 6   2rr\ 
cos  (-  +  y) 


where  cos  6 


=  -  J  -^L  ( 


kJ3     3Pf '  2 


3Pf 


=  0,  this  reduces  to 


=  -M. 


3H       '    TT 

Note  3:   For  definitions  of  coefficients,  see  Table  2.U.I. 
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Fig.  2.8.1  BASIC  ACTUATOR  PERFORMANCE  CURVES  FOR  LIFTING 
ROTOR  IN  Pt  D,  L  SYSTEM. 
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2 .9  The  Idealized  Rotating  Finite  Torque  Actuator 

The  ideal  actuator  mainly  considered  in  this  paper  is  essentially  a  zero 
torque  device.   In  section  9>  however,  a  detailed  analysis  is  presented  of  an 
actuator  which,  although  ideal  in  all  other  respects,  operates  at  finite  torque 
and  at  finite  shaft  speed.   Such  a  device  will  always  operate  somewhat  less 
effectively  than  the  ideal  zero  torque  actuator,  although  at  high  shaft  speed 
and  low  torque,  the  difference  will  be  small.   In  dealing  with  the  finite  torque 
actuator,  it  turns  out  to  be  advantageous  to  revert  to  the  use  of  the  conventional 
performance  coefficients  C  ,  C  and  J  based  on  p,  N  and  D  as  dimensional  reference 
parameters.   Moreover,  the  relations  among  C  ,  C  and  J  can  be  expressed  paramet- 
rically  in  terms  of  certain  slipstream  parameters  as  defined  in  Table  2.9-1. 

The  analysis  of  section  9  shows  that  the  torque  and  rotation  effects  lead 
to  the  creation  of  peripheral  velocities  in  the  slipstream  leaving  the  actuator. 
A  solution  is  found  for  these  peripheral  velocities  which  maintains  the  axial 
velocity  distribution  unchanged  from  what  it  would  be  for  the  ideal  zero  torque 
case.   In  other  words  the  axial  velocities  at  the  disc  and  in  the  remote  slip- 
stream remain  constant  and  uniform  as  before.  The  slipstream  is  found  to 
consist  of  two  distinct  regions,  namely,  a  small  inner  core  of  radius  r  in 
which  the  flow  remains  purely  axial,  and  an  outer  annulus  in  which  the  flow 
is  of  the  free  vortex  type. 

The  chief  results  of  this  analysis  are  summarized  in  Table  2.9.1-   Notice 
that  this  analysis  of  the  rotating  actuator  is  restricted  to  the  simple  modes 
P  and  W  only;  the  extension  of  the  analysis  to  include  all  modes  appears  to  be 
feasible  but  has  not  been  undertaken  in  this  paper. 

In  principle,  the  equations  of  Table  2.9.1  suffice  to  fix  a  family  of 
curves  of  C^  versus  J,  with  C_  as  the  parameter  whose  value  changes  from  one 
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Table  2.9.1  Conventional  Performance  Coefficients  of  Idealized 
Finite  Torque  Actuator  in  Terms  of  Slipstream  Parameters 

Dimensional  Slipstream  Parameters 

R  =  Slipstream  radius 

r  =  Slipstream  radial  co-ordinate 

(0  =  Shaft  angular  velocity  (rad/sec) 

V  =  Axial  velocity  of  undisturbed  flow 
Vc  =  Slipstream  axial  velocity  (constant) 

V  =  Slipstream  peripheral  velocity 

V  r  =  T  Slipstream  vortex  strength  (Constant  over  the  region 

rQ  *  r  £R.)  (2.9-1) 

Dimensionless  Slipstream  Parameters 

-=■     =  h     Slipstream  advance  ratio  (2.9-2) 

uiR 

(V  -  V) 

=  A  Slipstream  loading  parameter  (2.9-3) 

■p 

— -x  =  Y  Slipstream  swirl  parameter  (2.9.*0 

uiR 

x 

-=j-  =  Tl   Relative  radius  of  non-swirling  core  (2.9.5) 

Range  of  Loading  Parameter 

.  2  (>E-£\K   £  A  fi  +(7i  +  h2  -  ^  (2.9-6) 


Positive  values  of  A  correspond  to  mode  P,  negative  values  to  mode  W.  Other 
modes  are  not  considered  in  this  particular  analysis. 

(Table  continues  on  next  page) 
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Table  2.9.1         Cont'd 


Auxiliary  Functions 


-7: 


/I  -  2   (h  +  £)   A 
(h  +  |)  A 


(2.9-7) 


A,2    .2 


1  +  i   (h  +  A)    A  +  i   (H  +  ^      ^   + 


V  =   C   (h  +  |)    A 


(2.9-8) 


*-'** 


(2.9-9) 


Parametric  Solutions   for  Conventional  Performance  Coefficients 


In 


* 


Tb 


^ 


♦  i  A2   Tfc2 


(2.9.10) 


(2.9-11) 


J  =  m 


(2.9.12) 


Static  Conditions  and  Maximum  Loading  Limits 

Under  static  conditions  (V  =  0),  the  above  relations  simplify  as  indicated 
in  section  9*13 ?  with  results  as  shown  in  Fig.  2.9-1. 

Under  maximum  loading  conditions,  they  simplify  as  indicated  in  section 
9-1^,  with  results  as  shown  in  Fig.  2.9.2. 


(End  of  Table  2.9. l) 
3h 


curve  to  the  next.   The  calculations  required  for  this  purpose  are  straight- 
forward but  tedious.  At  the  time  of  writing,  they  still  remain  to  be  carried 
out.  When  completed,  the  resulting  family  of  curves  will  define  idealized 
performance  limits  which  provide  a  suitable  standard  of  comparison  for  the 
performance  of  any  real  single  rotation  actuator.  By  application  of  the 
rationale  explained  more  fully  in  section  9>  these  results  may  also  be  used 
to  evaluate  the  idealized  rotation  power  loss  P  for  any  prescribed  operating 
condition. 

Under  static  conditions  (V  =  0),  the  general  solution  summarized  in  Table 
2.9.1  reduces  to  the  special  form  given  in  section  9-13-   The  corresponding 
result  is  shown  in  Fig.  2.9-1  of  this  section  in  the  form  of  a  curve  of  the 
quanties  C_  and  C_/Cp   plotted  versus  Cp.  The  intercept  at  Cp  =  0  corresponds  to 
the  limiting  performance  at  zero  torque,  and  the  drop-off  of  this  relative 
thrust  parameter  with  increasing  power  coefficient  Cp  shows  the  effect  of 
rotation  losses. 

Another  special  case  of  theoretical  interest  is  that  corresponding  to  the 

maximum  loading  theoretically  attainable,  as  analyzed  in  section  9«1^.   The 

essential  results  are  summarized  here  in  Fig.  2.9.2  which  shows  C  ]    and 

VMAX 
Cp\    plotted  as  functions  of  J.   Note  that  this  absolute  performance  limit 

/MAX 
is  deduced  purely  from  fundamental  dynamical  principles  and  without  any  reference 

whatever  to  any  empirical  data  or  to  such  complex  and  highly  variable  details 

of  the  flow  field  as  local  lift  coefficients,  separation  effects  or  the  like! 

These  latter  effects  can  and  do  reduce  the  performance  of  any  actual  device 

well  below  the  limits  shown  in  Fig.  2.9-2.   The  essential  point,  however,  is 

that  there  is  no  feasible  way  to  raise  the  actual  performance  above  these 

limits . 
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It  is  also  noteworthy  that  Table  2.9.I  summarizes  the  performance  of  the 
ideal  finite  torque  actuator  in  a  comparatively  simple  algebraic  form.   In 
general  these  results  yield  information  that  goes  well  beyond  what  can  be 
deduced  from  the  usual  elementary  analysis  which  simply  neglects  rotation  effects 
altogether . 
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2.10  Stability  of  the  Ideal  Lifting  Rotor 

In  this  subsection,  we  summarize  the  dynamic  response  of  an  ideal 
lifting  rotor  (P  =  0)  to  small  perturbations  from  equilibrium  during  vertical 
ascent  or  descent.  At  the  equilibrium  state,  lift  L  is  equal  to  the  rotor 
craft  weight  W,  but  lift  is  subject  to  small  perturbations  6L  about  the  equilib- 
rium value . 

It  is  advantageous  in  this  application  to  choose  p,  D  and  W  as  dimensional 
reference  parameters.   The  corresponding  dimensionless  coefficients  of  power 
COTl,   climb  rate  C0^,  and  lift  C_  are  then  as  defined  in  Eqs .  (2.10.1)  through 

jr  Jn  jli 

(2.10.3)  of  Table  2.10.1.   It  also  simplifies  matters  to  express  time  in  terms 

of  the  dimensionless  coordinate  9  as  defined  in  Eq.  (2.10.U). 

The  inherent  stability  characteristics  of  the  lifting  rotor  are  expressed 

/ac  \    /ac  \ 

by  the  partial  derivatives   -_    and    -*   which  must  be  evaluated  along  the 

V^SH/    \3C3L/ 
equilibrium  operating  line  C_  =  1.  For  simplicity,  these  derivatives  are 

denoted  as  X  and  Y,  respectively,  and  are  known  functions  of  C  •,  as  summa- 

rized  in  Eqs.  (2.10.5)  through  (2.10.8). 

The  power  input  to  the  rotor  or,  in  descent,  the  power  output  from  the 
rotor  may  either  be  held  constant,  or  may  be  modulated  by  a  control  system 
in  response  to  sensed  values  of  the  perturbation  acceleration  and/or  perturbation 
velocity,  according  to  Eq.  (2. 10. 9).  Constants  p,  and  v  express  these  two 
components  of  the  power  response. 

The  final  dynamic  response  of  the  rotor  craft  to  a  small  perturbation 
is  then  exponential  in  character  as  shown  in  Eq.  (2.10.10).  A  stable  response 
corresponds  to  a  negative  value  of  the  exponential  factor  (3. 

The  solution  for  (3  is  given  by  Eq.  (2.10.11).   This  is  seen  to  depend 

on  the  functions  X  and  Y  which  describe  the  inherent  stability  characteristics 

of  the  rotor  itself,  and  upon  the  constants  u.  and  v  which  characterize  the  power 

control  system,  if  any. 
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Table  2.10.1  Stability  of  Idealized  Lifting  Rotor 
in  Vertical  Ascent  or  Descent 


(Pf  =  o) 


Dimensionless  Coefficients 


'3P 


-  (2.10.1)  C3L=|  (2.10.3) 


1 

D* 


—2—   (2.10.2)   9  =  ±—   (2.10.1+) 

Aj  1   /w 


1  A[ 
D  J   P 


gP  V  P 


Actuator  Derivatives  at  Equilibrium   (C_     =  l) 

Jl; 

Modes   P,  W  and  R 


x-^P  =  i  + 
3C3H 


H 


Y_^3P_^H  + 
BGT     "2 


'H\       ,    2T 


V\2 

2  /        rr  cp 


ffi/T     +^ 


(2.10.5) 


(2.10.6) 


Mode  T 


x  .  Zk  I  -±  T  J  1 1 


dC 


3H 


^ 


H 


(2.10.7) 


BC3L  CH 


(2.10.8) 


(Table  continues  on  next  page) 
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Table  2.10.1    Cont'd 


Control  System  Coefficients 


6C3P  =  "  *  ^   'L5C3H 


5C 


3H 


(2.10.9) 


m  =  perturbation  acceleration  response  constant  (if  any)  ^  0 
v  =  perturbation  velocity  response  constant  (if  any)  ^  0 

Dynamic  Response  to  a  Perturbation 


*3H    "  |  6C3h 


.09 


J    L     J 


where 


(2.10.10) 


0  =  - 


X  +  v 

Y  +  p, 


(2.10.11) 


Algebraic  Signs  of  Actuator  Derivatives 


X  >  0  for  all  values  of  COTv 

Y>0    ifC3ft>-j? 

Y  <  0   if  C  •  <  -  -2- 
3H     ^ 


Stability 

Range  of  stable  response  for  p  =  v  =  0 


(2.10.12) 
(2.10.13) 


0  <  0  for   -  -4s  <  C~v  <  +  * 

VTT      J"11 


(2.10.1U) 


Stability  Augmentation  - 

1.  Increasing  p  increase  stable  range 

2.  Increasing  v  - 

a)  does  not  affect  stable  range 

b)  increases  response  rate  (whether  stable  or  unstable) 

(End  of  Table) 
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Calculations  reveal  that  the  rotor  derivatives  X  and  Y  are  as  shown  in 
Fig.  2.10.1.   This  leads  to  the  conclusion  that  with  |i  =  v  =  0,  the  system 
is  inherently  stable  over  the  range  shown  in  Eq.  (2.10.1U).   Inspection  of 
Eq.  (2.10.11)  then  shows  that  the  range  of  stable  operation  can  be  extended 
at  will  by  increasing  the  acceleration  control  constant  |i.   Increasing  the 
velocity  control  constant  v  does  not  affect  the  stability  range  itself,  but 
augments  the  response  rate,  whether  stable  or  unstable. 

The  exponential  response  rate  g  which  fixes  the  degree  of  stability  or 
instability  is  plotted  versus  climb  rate  C  •  for  several  cases  of  interest  in 
Fig.  2.10.2. 

It  is  to  be  expected  that  the  stability  characteristics  of  any  real  lifting 
rotor  will  generally  be  considerably  less  favorable  than  the  stability  of  the 
ideal  rotor  as  summarized  here. 
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2.11  Auxiliary  Figures  of  Merit 

The  purpose  of  this  brief  subsection  is  simply  to  call  attention  to  a 
number  of  auxiliary  quantities  which  do  not  occur  in  the  main  theoretical 
development  but  which  nevertheless  provide  useful  figures  of  merit  in  certain 
applications.   These  include  the  concepts  of  efficiency,  power  effectiveness 
and  force  effectiveness  as  defined  in  section  11. 
3-  Basic  Assumptions 

The  purpose  of  this  analysis  is  not  to  predict  as  accurately  as  possible 
the  performance  of  any  particular  real  actuator  under  specified  conditions. 
It  is  rather  to  establish  as  simple  as  possible  a  theoretical  performance  limit 
which  any  real  actuator  might  approach  but  which  none  could  ever  exceed.   In 
general,  the  performance  of  any  actual  device  involves  physical  mechanisms  so 
complex  that  their  effects  cannot  be  adequately  assessed  by  theoretical  methods 
alone;  it  is  necessary  to  refer  also  to  the  results  of  experiment.  Fortunately, 
the  determination  of  a  simple  performance  limit  is  a  much  easier  task.  For 
this  purpose  it  is  not  only  permissible  but  mandatory  to  neglect  all  complex 
but  secondary  effects,  and  to  concentrate  attention  exclusively  upon  the 
irreducible  essentials  of  the  phenemena  under  study. 

For  this  reason  we  treat  the  fluid  as  incompressible  and  inviscid,  we 
treat  the  flow  field  as  possessing  polar  symmetry  about  the  axis  and  as  being 
steady,  and  we  neglect  possible  hydrodynamic  instability  and  turbulence  effects. 
As  a  result  of  passing  through  the  actuator,  the  fluid  ultimately  undergoes  a 
change  of  axial  velocity  A  V.  We  assume  that  this  change  develops  asymptotically 
within  the  slipstream  as  the  fluid  moves  away  from  the  actuator  and  that,  once 
established,  the  slipstream  persists  indefinitely.   In  particular,  for  any  of 
the  flow  patterns  involving  flow  reversal,  or  involving  co-axial  streams 
moving  at  different  velocities  and  possibly  even  in  opposite  senses,  we  assume 
that  these  reversed  or  co-axial  streams  persist  indefinitely  without  breaking  down. 
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For  any  real  actuator  operating  in  a  real  fluid,  of  course,  the  slip  stream 
eventually  disappears  owing  to  the  effects  of  hydrodynamic  instabilities, 
turbulence  and  viscous  dissipation.   These  effects  are  particular  important 
for  cases  involving  flow  reversal  or  co-axial  streams,  especially  if  the  streams 
move  in  opposite  senses.  Nevertheless,  even  for  real  flows,  it  is  usually 
possible  to  locate  upstream  and  downstream  cross-sections  which  are  sufficiently 
close  to  the  actuator  that  the  slipstream  is  still  well  defined  and  yet  far 
enough  away  from  the  actuator  that  the  pressure  in  the  stream  is  substantially 
equal  to  ambient  pressure.  For  all  practical  purposes,  conditions  at  such  cross- 
sections  are  comparable  to  conditions  "at  infinity"  in  the  idealized  flow  pattern. 
Hence  the  neglect  of  slipstream  dissipation  effects  is  not  as  drastic  an  ideal- 
ization as  it  might  seem  at  first  glance.   Moreover,  the  usefulness  of  this 
approach  has  long  been  established  in  connection  with  the  simple  modes  P  and 
W.   In  this  paper  we  are  extending  the  same  method  to  the  other  modes  R  and  T. 
Admittedly,  these  latter  modes,  especially  mode  T,  are  more  complex;  nevertheless, 
for  the  special  purpose  of  establishing  a  simple  performance  limit,  the  method 
appears  to  be  fully  justified. 

The  idealized  analysis  also  assumes  that  the  velocities  are  purely  axial  at 
all  reference  cross-sections.  Of  course,  any  real  propeller  or  other  rotating 
actuator  always  exerts  some  torque  upon  the  fluid  and  therefore  sets  up  corre- 
sponding peripheral  velocity  components  in  the  slipstream.  Creation  of  the 
kinetic  energy  associated  with  these  components  absorbs  shaft  power.   The 
resulting  reduction  of  pressure  over  much  of  the  slipstream  also  adversely 
affects  the  net  thrust.  Hence  the  peripheral  velocity  components  degrade 
performance  and  are  therefore  eliminated  from  the  idealized  model,  except  of 
course  in  section  9  where  this  matter  is  considered  in  detail. 
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Similarly,  a  real  propeller  or  non-optimal  actuator  may  create  non-uniform 
axial  velocities  in  the  far  slipstream.  Thrust  is  proportional  to  the  mean  axial 
velocity.   Power  required  is  proportional  to  the  mean  square  of  the  axial  velocity. 
Clearly,  the  ratio  of  useful  thrust  delivered  to  power  required  will  be  a  max- 
imum only  if  the  axial  velocity  is  uniform  across  the  slipstream.  Consequently, 
the  analysis  of  the  idealized  actuator  assumes  uniform  axial  velocities  across 
the  far  slipstream  and  across  all  reference  stations.   Pressures  are  also 
taken  as  uniform  across  these  stations. 

In  the  transition  mode  T  which  involves  two  oppositely  moving  co-axial 
streams,  the  velocities  and  pressures  are  taken  as  uniform  across  each  stream 
separately  at  each  station;  however,  the  values  of  velocity  and  pressure  may 
in  general  be  different  for  the  two  separate  streams  at  a  given  cross -section. 

In  the  classical  momentum  analysis  of  the  simple  modes  P  and  W,  it  is 
assumed  implicitly  that  the  net  axial  force  between  the  fluid  which  passes 
through  the  actuator  and  the  ambient  flow  is  zero  or  at  least  that  it  is  negli- 
gible compared  with  the  axial  force  exerted  at  the  actuator  disc  itself.   The 
corresponding  assumption  is  here  extended  to  modes  R  and  T.   In  other  words, 
the  net  axial  force  between  any  two  co-axial  streams  is  assumed  to  be  negligible 
in  comparison  with  the  axial  force  exerted  on  either  stream  at  the  actuator 
disc  itself.   The  only  exception  to  this  occurs  in  section  9  which  deals  with 
the  finite  torque  actuator;  in  this  case  it  is  necessary  to  admit  a  net  axial 
force  across  the  interface  between  the  nonrotating  core  and  the  free  vortex 
flow  region. 
k.     Additional  Nomenclature  for  Main  Text 

All  symbols  used  in  section  2  have  previously  been  defined  in  subsection  2.1. 
All  other  symbols  occurring  in  this  paper  and  not  previously  defined  in  subsection 
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2.1  are  summarized  below.  Again  the  definitions  are  arranged  according  to 

the  section  and  in  the  approximate  order  in  which  they  are  introduced  in  the 

text. 

Section  3 

AV  =  (V„    -   V)  =  net  overall  increase  in  axial  velocity  occurring  across 

the  actuator,  positive  if  in  the  same  sense  as  the  velocity  V  of  the 

undisturbed  fluid,   ft /sec 

Note:  For  the  important  special  case  of  static  thrust  defined  by  the  condition 

V  =  0,  the  velocities  V,,  V  and  AV,  and  the  force  F  exerted  by  the  actuator 

upon  the  fluid  are  all  in  the  same  sense  and  all  have  the  same  sign  which  is 

normally  taken  as  positive. 

Section  5 

0,  u,  d,  s  =  symbols  or  subscripts  denoting  the  following  locations, 

respectively:   far  upstream  side  of  actuator  disc,  downstream 

side  of  actuator  disc,  and  slipstream  far  downstream  of  actuator 

V,  V  ,  V, ,  V„  =  respective  axial  velocities  at  the  above  locations,  with 
u   d   S 

V  =  V..   ft/sec 
u    d 

p,  p  ,  p  ,  p  =  respective  pressures  at  the  above  locations,  with  p  =  p. 

lbf/ft2 
E  =  total  energy  per  unit  mass  of  fluid  upstream  of  actuator,  ft  lbf/slug 
AE  =  change  in  total  energy  per  unit  mass  of  fluid  produced  by  actuator. 

ft  lbf/slug 
m  =  mass  flow  rate  through  actuator  disc,   slugs /sec 

Pa  =  power  lost  to  slipstream  by  ideal  zero  torque  actuator,   ft  lbf/sec 
Section  6 

C    =  ideal  power  coefficient  in  p,  D,  V  system.  -- 


\  =  —  =   velocity  change  ratio, 
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U  =  symbol  designating  that  branch  of  the  formal  solution  which  represents 

physically  unstable  operation. 

S,  0,  M,  N,  U,  V,  A,  B,  C,  S'  =  symbols  designating  significant  points  on 

ideal  operating  envelope  as  shown  in  Fig.  6.2 

Section  7 

x  =  fraction  of  actuator  disc  area  operating  in  the  windmill  mode  during 

operation  in  the  transition  range. 

/dClp\ 

S  =  I  — )  =  slope  of  ideal  operating  curve 

VC1F/ 

A,  C  =  subscripts  referring  to  operation  in  windmill  and  reversed  modes, 

respectively,  at  points  A  and  C  of  envelope. 

Section  8 

V  =  ultimate  axial  velocity  of  propagation  of  vorticity  generated  at 

actuator  disc,   ft/sec 

Section  9 

Subsection  9»2 

P.  '  =  shaft  power  of  ideal  finite  torque  actuator,   ft  lbf/sec 

P   =  idealized  rotation  power  loss,  always  non-negative,   ft  lbf/sec 

P   =  residual  power  loss,  always  non-negative,  ft  lbf/sec 

Subsection  ^.h 

p  =  pressure  at  edge  of  slipstream  at  radius  R  (equal  to  pressure  of 

p 
undisturbed  flow),   lbf/ft 

p„  =  pressure  within  slipstream  at  radius  r  <  R  (lower  than  pressure  of 

undisturbed  flow),  lbf/ft 

V  =  peripheral  velocity  in  slipstream,  ft/sec 
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Subsection  9»5 

o 
pqn  =  pressure  in  central  non-swirling  core  region  of  slipstream.  Ibf/ft 

x 
7|  =  —  =  slipstream  relative  radius. 

Subsection  9«1Q 

A^.  ,  Cma-y>  Vmay  =  values  of  the  quantities  A,  £,  Y  under  conditions  of 

maximum  loading  attainable  in  an  idealized  rotating 

actuator. 

Subsection  9*15 

B,  =  relative  angle  of  flow  at  radius  r,  on  downstream  side  of  actuator 
Md  d 

as  defined  in  Eqs.  (9.15.I)  or  (9.I5.U).  rad 
3  =  relative  angle  of  flow  at  relative  radius  T]  in  remote  slipstream 
as  defined  by  Eq.  (9.15-5).   rad 

TL  =  mean  radius  ratio  as  defined  by  Eq.  (9. 15. 6). 

Section  11 

T)  =  overall  efficiency  of  real  actuator  as  defined  by  Eq.  (ll.l.). 

T].  =  efficiency  of  ideal  actuator  as  defined  by  Eq.  (11.1.2). 

e  =  rotation  power  effectiveness  as  defined  by  Eq.  (11.2.1). 

e  =  residual  power  effectiveness  as  defined  by  Eq.  (11.2.2). 

e  =  overall  power  effectiveness  as  defined  by  Eqs.  (11.2.3)  or  (11. 5). 

e^  =  force  effectiveness  as  defined  by  Eq.  (11.6.2). 
r 

e_,  =  braking  effectiveness  as  defined  by  Eq.  (11. 6. 3). 

D 

Cn_\    =  maximum  value  of  C._  for  an  ideal  actuator  at  a  specified  value 
)  MAX 

ofClp.   - 

Cn„\    =  minimum  value  of  Cni_,  for  an  ideal  actuator  at  a  specified  value 

1fJmax  1F 

of  Clp.   — 
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5-   Momentum  and  Energy  Analysis  of  Idealized  Zero  Torque  Actuator  in  Modes 

P,  W  and  R 

For  all  three  of  the  flow  configurations  P,  W  and  R  considered  in  this 
section,  there  are  four  major  reference  stations  within  each  stream,  designated 
respectively  by  subscripts  0,  u,  d  and  s.   The  velocities,  pressures  and  total 
flow  energies  at  these  stations  are  designated  as  shown  in  Table  5.1. 


Table  5.1  Velocity,  Pressure  and  Energy  at  Principal  Stations 


Total 

Energy 

Static 

Per  Unit 

Station 

Location 
Far  upstream 

Velocity 
V 

Pressure 
P 

Mass 

0 

E 

u 

At  actuator, 
upstream  side 

Vu  =  Vd 

Pu 

E 

d 

At  actuator, 

Vd 

Pd 

E  +  IE 

downstream  side 

s 

Slipstream 

Vg  =  V  +  AV 

P 

E  +  t$ 

It  should  be  noted  that  the  terms  upstream  and  downstream  in  Table  5.1 
and  elsewhere  in  this  discussion  always  refer  to  the  effective  flow  through 
the  actuator,  not  to  the  ambient  flow  field.   Thus  for  reversed  flow  in  mode 
R,  the  upstream  station  u  is  on  the  aft  side  of  the  actuator  disc!  Also  note 
that  by  reason  of  continuity,  the  axial  velocity  must  be  equal  at  stations  u 
and  d  and  that  the  single  symbol  V,  is  therefore  used  for  both. 

Bernoulli's  equation  can  be  applied  firstly  between  stations  0  and  u  and 

secondly  between  stations  s  and  d.  Thus  for  any  single  mode  we  thereby  obtain 

2 

T72      P       V/ 

E=P+f  =  -2i  +  -d-  (5.1.1) 

p    2     p     2 
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e  +  m  -  £  +  (v  Vv)g  -^*\-  (5.i.2) 

p       2        p    2 

Upon  subtracting  ( 5.1.1)  from  (5.1.2)  and  rearranging,  we  obtain  the  difference 
in  total  energy  and  in  static  pressure  across  the  actuator  disc  in  the  form  - 

p/E  =  (pd  "  Pu)  =  P  (v  +  -f)  AV  (5.2) 

The  present  analysis  is  restricted  to  unshrouded  configurations  only. 
Thus  the  whole  thrust  F  is  produced  by  the  above  static  pressure  difference 
acting  on  the  actuator  disc  area.  However,  the  sense  of  the  thrust  depends 
on  whether  the  sense  of  the  flow  through  the  disc  is  positive  or  negative. 

We  therefore  express  the  result  in  the  generalized  form 

p  =  ^  pd2  (v  +  f )  m  (5.3) 

where  the  operator  cp  assures  that  the  force  always  has  the  correct  algebraic 
sign. 

A  second  independent  expression  for  the  thrust  can  also  be  found  by 
applying  the  momentum  theorem  between  stations  0  and  s .  For  this  purpose  we 
need  the  absolute  mass  flow  rate  through  the  actuator  which  can  be  expressed 
in  the  form 

Consequently  the  required  thrust  may  be  finally  written  as 

F=^pD2VdAV  (5-5) 

Comparison  of  Eqs .    (5-3)  and   (5.5)  now  reveals  that  algebraically 

Vd  =  V  +  f  (5.6) 
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For  the  idealized  actuator,  the  shaft  power  P.  goes  entirely  into 

changing  the  total  energy  of  the  fluid  which  passes  through  it.  Thus  the 
algebraic  shaft  power  P.  equals  the  product  of  the  algebraic  energy  change 

M   as  expressed  in  Eq.  (5.2)  with  the  absolute  mass  flow  rate  m  as  given  by 

Eq.  (5-^).  Therefore,  using  Eq.  (5-6)  to  eliminate  V,,  we  finally  obtain 
for  the  power 

2 


P. 


i=?P^  V  +  f   AV  (5.7) 


At  this  point  the  key  relations  applicable  to  modes  P,  W  and  R  can  be 
summarized  as  follows : 


vs 

— 

V  + 

AV 

vd 

= 

V  + 

AV 
2 

CP 

= 

(5.8.1) 
(5.8.2) 


and 


F=?PD2  (v  +  f)  AV  (5.9.1) 

P±  =  9  p  D2  (v  +  m      AV  (5.9-2) 

Eqs.  (5.9)  define  the  fundamental  relation  among  the  principal  quantities 
p,  A,  V,  F  and  P.  in  terms  of  the  auxiliary  parameter  AV.  While  it  is  not 
mandatory  to  do  so,  it  is  in  some  respects  convenient  to  eliminate  this 
secondary  parameter  AV  between  the  two  relations  (5.9).   The  result  may  be 
written 

F3  =  m      D2  P>  (P  .  FV)  (5.10.1) 

2  K     1   1 
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We  note  that  this  equation  is  linear  in  V,  quadratic  in  P.  and  cubic  in 

F.  Hence  explicit  solutions  for  V  or  P.  may  be  readily  obtained,  but  an 

explicit  solution  for  F,  while  entirely  feasible,  turns  out  to  be  more  involved. 

We  shall  not  undertake  the  cubic  solution  for  F  at  this  point.  However,  solving 

Eq.  (5.10.1)  for  V  gives 

P..    n      /  „2   \ 

(5.10.2) 


Likewise,  solving  Eq.  (5.10.1)  for  P.  we  obtain  two  roots.  Nevertheless,  by 
use  of  the  operator  a.,  the  solution  may  be  expressed  unambiguously  as  follows. 

It  is  shown  in  a  later  section  that  the  formal  analytical  solution  defined 
by  Eqs.  (5.9) >  or  by  Eqs.  (5.10.1)  or  (5.10.2),  consists  of  four  distinct" branches 
Three  of  these  correspond  respectively  to  modes  P,  W  and  R.  The  fourth  branch 
is  extraneous;  it  represents  a  theoretical  solution,  mode  U,  which  is  actually 
unstable  for  any  unshrouded  configuration.  One  advantage  of  the  particular 
form  of  solution  defined  by  Eq.  (5. 10. 3)  is  that  by  reason  of  the  selective 
character  of  the  operator  a.  ,  this  extraneous  branch  is  eliminated. 

It  is  also  instructive  to  analyze  the  basic  momentum  energy  relations  with 
respect  to  axes  which,  instead  of  being  fixed  with  respect  to  the  actuator 
itself  as  assumed  heretofor,  are  fixed  instead  in  the  remote  fluid.   In  this 
axis  system,  the  actuator  is  moving  with  velocity  V.  The  reaction  force  of 
amount  F  exerted  by  the  fluid  upon  the  actuator  therefore  delivers  a  thrust 
power  of  amount  FV.  Of  course,  this  may  be  positive  or  negative  thrust  power, 
depending  on  the  sign  of  F. 
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In  addition,  the  velocities  at  stations  0  and  s  with  respect  to  the  new 

axes  are  now  zero  and  AV,  respectively.   The  corresponding  total  energies  per 

,2" 


unit  mass  are  *-  and 

P 


2  +  iMl 

P     2 


,  respectively.  Hence  the  net  change  in  total 

2 


energy  per  unit  mass  is  the  difference  "between  these  two  quantities  or  ' 


2 
This  is  simply  the  kinetic  energy  per  unit  mass  of  the  slipstream  with  respect 

to  the  new  axes.  Since  the  slipstream  flow  is  being  "generated"  at  the  rate 

of  m  units  of  mass  per  unit  time,  the  creation  of  this  kinetic  energy  in  the 

2 

slipstream  consumes  power  at  the  rate  P  =  m  '  0   .  Of  course  the  power  Pa 
imparted  to  the  slipstream  represents  a  loss  of  useful  power  which  is  basic 
and  unavoidable,  even  for  an  actuator  which  is,  by  definition,  completely 
idealized.   Moreover,  since  the  thrust  power  FV,  the  power  Pq  lost  to  the 
slipstream,  and  the  shaft  power  P.  must  satisfy  an  overall  energy  balance,  we 
may  write  at  once  that 

PS  =  (Pi  "  FV)  =  A  "^P~"  *  °  (5.11.1) 

By  utilizing  Eqs.    (5.*+)   and   (5-5)?   this  relation  can  readily  be  amplified  to 
reveal  that 

Pg  .    L     .  Fvj    ,lf -.A  i^tz0  (5.11.2) 

F  AV 
This  relation  in  terms  of  the  quantity  -5—  can  also  be  obtained  independently 

by  dividing  Eq.  (5. 9-2)  by  (5-9.1)  and  rearranging.   This  agreement  provides 

further  confirmation  of  Eqs.  ( 5.11.1)  and  (5.11.2).   It  also  shows  that  the 

power  analyses  in  the  two  different  axis  systems  are  mutually  consistent. 

It  is  significant  that  the  separate  quantities  P.  and  FV  can  be  either 

positive  or  negative,  but  that  the  difference  (P.  -  FV)  must  always  remain 

non-negative.  As  explained  above,  this  quantity  represents  an  irreducible 

minimum  loss  of  power  which  is  unavoidable  even  with  an  ideal  actuator. 
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Naturally,  for  any  real  actuator  operating  with  corresponding  values  of  p,  D,  V 
and  F,  the  actual  power  loss  must  be  even  greater  than  this  basic  minimum. 

By  rearranging  Eqs.  ( 5.10.1)  and  (5. 10. 3)  we  can  also  find  two  alternative 
and  equivalent  expressions  for  the  basic  power  loss  in  terms  that  are  frequently 
more  convenient  than  the  original  form  (5.11.2).  They  are,  respectively, 


>P  D  P. 


and 


h-m-U  m  +1:4  -w«         «.a.» 


Once  again,  the  sign  of  the  radical  in  Eq.  (5.11.U)  is  specified  by  the 
operator  a.  in  such  a  way  as  to  eliminate  the  physically  extraneous  branch  of 
the  solution. 

Notice  that  the  above  results ,  Eqs.  ( 5-8.1)  through  (5-ll.l+)  have  been 
established  in  a  form  which  does  not  involve  division  by  V  in  any  term.  Hence 
the  solution  for  the  important  special  case  of  static  thrust  may  be  found 
simply  by  substituting  V  =  0  into  these  equations.   In  this  case  we  also  find 
that  since  thrust  power  vanishes, 

\   =  Ps  (5.12) 

It  may  be  remarked  in  passing  that  a  fan,  unlike  a  propeller,  windmill 
or  lifting  rotor,  is  utilized  for  the  very  purpose  of  creating  a  moving  slip- 
stream. Hence  the  power  Pq  imparted  to  the  slipstream,  which  must  be  regarded 
as  a  dead  loss  for  these  other  three  devices,  becomes  the  principal  useful 
output  for  the  fan!  Apart  from  this  single  change  in  purpose  and  viewpoint, 
the  analysis  of  the  idealized  fan  is  exactly  the  same  as  that  given  here  for 
these  other  actuators. 

56 


The  fundamental  result  of  this  section  is  the  basic  actuator  relation, 
Eq.  (5.IO.3).  Alternatively,  we  may  choose  either  of  the  equivalent  forms 
(5.10.1)  or  (5.10.2).  The  same  basic  relation  is  also  defined  by  the  pair 
of  parametric  equations,  (5.9«l)  and  (5.9-2),  in  which  AV  plays  the  role  of 
auxiliary  parameter. 

On  the  other  hand  for  the  exceptional  case  of  a  fan,  the  basic  relation 
between  output  Pa  and  input  P.  is  summarized  parametric ally  by  the  pair  of 
equations,  (5.11.3)  and  (5.IO.3),  in  which  F  now  plays  the  role  of  auxiliary 
parameter. 

The  foregoing  fundamental  results ,  in  any  or  all  of  the  above  forms ,  may 
be  further  simplified  by  non-dimensionalization.   This  process  reduces  the 
number  of  parameters  finally  involved  by  three.   The  results  previously  sum- 
marized in  section  2.0  were  developed  in  this  way. 

Eq.  (5.10.1)  which  is  cubic  in  F,  is  not  solved  explicitly  for  F  in  this 
section.  However,  the  cubic  solution  has  been  worked  out  in  detail  and  pre- 
sented in  connection  with  several  of  the  non-dimensional  forms  summarized  in 
section  2.0. 
6.  Stable  and  Unstable  Branches  of  the  Formal  Solution 

It  is  useful  at  this  point  to  non-dimensionalize  the  foregoing  results  in 
terms  of  the  dimensional  reference  parameters  p,  D,  V.  For  this  purpose  we 
define  the  following  dimensionless  coefficients,  namely, 


'IF 


D2V* 


(6.1.1) 


P. 

Cip,  =  i-5T  (6.1.2) 

1Pl    p  D2  V3 


X=f  (6.1.3) 
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In  terms  of  these  coefficients,  the  detailed  formal  solution  represented 
by  Eqs.  (5. 8.1)  through  (5.11.U)  may  be  summarized  in  dimensionless  form  as 
follows . 


'ir)  =  1  +  2  <6-2-2> 


¥  =  -) ¥;  (6.2-3) 


°1F  "  ¥  I1  +  iU  <6-3-i> 


-&f*i2* 


C1Pr  -  if  (1  +  #1     X  (6.3.2) 


cif3  ■  f  cm  (cm  "  °!,)  <*>M> 


^-M^yw-i^  M.« 


The  last  three  of  these  equations  as  given  below  represent  alternative 
expressions  for  the  dimensionless  power  loss  to  the  slipstream. 


(6.5.1) 


"    flPi  "  cif] 


2 

2  TTCp  Clpi 
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i  J\2J       +  tkpC1F   ■  T" 


(6.5.3) 
(\\         /v  A 

The  velocity  ratios  \— j   and  (— )  as  expressed  by  Eqs.  (6.2)  are  plotted 
as  functions  of  parameter  \   in  Fig.  6.1.  The  corresponding  relation  between 
C^  and  Clpi,  as  defined  parametric  ally  by  Eqs.  (6.3)  or  explicitly  by  Eq. 
(6.U.1),  is  plotted  in  Fig.  6.2. 

Study  of  these  charts  discloses  the  existence  of  four  distinct  branches  of 
the  formal  solution  which  are  correlated  with  the  previously  defined  modes  as 
follows : 

Branch  SO   -  Mode  P  (Propulsion) 
Branch  OMN  -  Mode  W  (Windmill) 
Branch  NUV  -  Mode  U  (Unstable) 


Branch  VCS'  -  Mode  R  (Reversed  Flow) 

To  see  that  branch  NUV  represents  a  physically  unstable  mode,  we  take  note 
of  the  following  facts.  For  this  branch,  parameter  \   lies  between  -1  and  -2. 
Eqs.  (6.2)  and  Fig.  6.1  then  show  that  the  flow  through  the  actuator  is  in  the 
same  sense  as  the  ambient  flow  while  the  final  flow  in  the  remote  slipstream 
is  reversed  with  respect  to  the  ambient  flow.   Thus  the  flow  reversal  takes 
place  after  the  fluid  passes  through  the  actuator.  Moreover,  Eqs.  (6.3)  and 
Fig.  6.2  show  that  C._  and  C-_.  are  both  negative  in  this  range.   In  other 

.Lr       -LitI 

words  the  actuator  sustains  a  drag  force  and  extracts  useful  power  from  the 
stream  which  subsequently  reverses  and  departs  in  a  sense  opposite  to  that 
of  the  ambient  flow. 

While  the  above  behavior  would  theoretically  be  possible  within  the  confines 
of  a  properly  designed  channel  with  physical  restraining  walls,  it  cannot  be 
sustained  in  the  unrestrained  environment  of  an  unshrouded  actuator.  The 
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mechanism  of  vorticity  transport  which  causes  this  flow  pattern  to  break  down 
is  discussed  in  some  detail  in  section  8.   It  suffices  here  merely  to  point  out 
that  a  flow  can  be  completely  reversed  with  an  unshrouded  actuator  only  by  the 
expenditure  of  shaft  power.  Actually,  the  flow  must  be  reversed  before  it  passes 
through  the  actuator,  and  must  therefore  enter  it  from  the  aft  side.  This 
requires  the  creation  of  a  sufficiently  intense  suction  on  the  aft  side  to  turn 
the  ambient  flow  inward  and  around  toward  the  disc.  Since  branch  NUV  does  not 
satisfy  these  various  requirements,  it  must  be  discarded  as  physically  irrelevent 
for  any  unshrouded  configuration. 

The  reader  can  verify  that  the  analytical  solution  in  the  specific  form  of 
Eq.  (6.U.2)  automatically  eliminates  the  above  extraneous  solution  through  the 
selective  action  of  the  operator  a. • 

We  note  in  passing  that  the  detailed  solution  which  is  non-dimensionalized 
in  this  section  in  terms  of  parameters  p,  D,  V  can  also  be  non-dimensionalized 
in  terms  of  other  reference  parameters,  for  example,  p,  D,  |pj  or,  in  the  case 
of  a  lifting  rotor,  p,  D,  L.  While  the  results  so  obtained  differ  in  algebraic 
form,  they  express  relationships  that  are  theoretically  equivalent  to  those 
summarized  above.  The  form  based  on  p,  D,  V  happens  to  be  most  convenient  for 
the  present  phase  of  the  analysis. 
7.  A  New  Solution  for  Transition  Range,  Mode  T 

Since  the  transition  from  mode  W  to  mode  R  cannot  actually  take  place 
along  the  unstable  branch  MJV  of  the  original  formal  solution,  we  are  compelled 
to  formulate  an  alternative  hypothesis  to  apply  over  this  range.  We  now  propose 

the  following  novel  yet  simple  solution  for  this  case. 

rr  2 
Imagine  the  total  disc  area  j-  D  to  be  subdivided  into  two  concentric  regions , 

p  TT    2 

namely,  an  inner  circle  of  area  x  £  D  and  an  outer  annulus  of  area  (l  -  x)  j-  D 
where  the  transition  parameter  x  lies  somewhere  in  the  range  between  0  and  1. 
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We  assume  that  the  flow  over  the  inner  region  corresponds  to  some  as  yet 
unspecified  operating  point  A  of  mode  W,  and  that  flow  over  the  outer  region 
corresponds  to  some  as  yet  unspecified  operating  mode  point  C  of  mode  R. 

Consider  the  following  relations. 
For  the  inner  flow  - 

-1  £  XA  £0  (7.1.1) 

/    X  \ 
C,«  -+"?(l  +  4'k  (7.1.2) 


2  )KA 

TT  /,      A 


cip«  =  +  ir  1  +  -    \  (7-1-3) 


For  the  outer  flow  - 


-  »  £  ^  £  -  2  (7.2.1) 


°ifc  "  "  f  f1  +  T-Vo  <7-2-2) 


cipic  "  -  f  f1  +  ¥)    xo  <7-2-3) 


The  corresponding  net  overall  values  for  the  combined  flow  are  found  simply 
by  adding  the  respective  contributions  from  the  inner  and  outer  streams.  The 
results  are 


C1F  "  X  C1FA  +  (1  "  X)  C1FC  (7'3-l) 


C1P1  =  X  ClPiA  +    (1   "  X)   ClPiC  (7'3'2) 
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Upon  eliminating  x  between  these  last  two  relations  we  obtain 


ClPi  ~  ClPiC  + 


n    -  p 
lPiA    lPiC 


1FA   ^1FC 


C1F  "  cifc 


(7.U) 


The  result  shows  that  for  fixed  values  of  \.  and  X^,  the  relation  between 

C,  ~.  and  C__  is  linear  over  the  transition  range. 
IPi      IF 

The  problem  that  remains  is  to  find  the  values  of  the  parameters  \.  and 
\  such  as  to  optimize  actuator  performance  over  the  transition  range.  This 
means  that  for  any  given  value  of  C-.„,   the  power  coefficient  C   ,  whether  positive 

Xr  JLr 

or  negative,  should  be  algebraically  as  small  as  possible.   Inspection  of  Fig. 
6.2  shows  that  this  requirement  is  satisfied  simply  by  choosing  as  the  transition 
line  the  straight  line  segment  AC  which  is  tangent  to  the  original  analytical 
curve  at  the  two  points  A  and  C  as  shown. 

In  this  connection  we  note  that  the  slope  of  the  original  curve  at  these 
two  points  can  be  expressed  in  the  form 


'd  C 


lPi- 


1  + 


d  C 


IF 


=  SA  = 


M 

2  / 


3X« 


1  + 


A 


TTT-^T 


(7.5-1) 


a^)(x>^ 


(1  +  \J 


(7.5.2) 


Furthermore,  the  slope  of  the  straight  line  segment  joining  points  A  and 
C  can  be  reduced  to  the  form 


S  = 


(7-5.3) 


6*+ 


Now  the  optimizing  condition  requires  simply  that 

s  =  sA  =  sc  (7.5.10 

Simultaneous  solution  of  Eqs .  (7-5)  gives  the  following  important  results  - 

S  =  SA  =  Sc  =  -1  (7.6.1) 


XA  =  -2  (^-f1)   =  -  0.9*5  (7.6.2) 


Xc  =  -2  f^rM  =  "  3.155  (7.6.3) 

Therefore  Eqs.    (7.1)   and   (7.2)   give  at  point  A  - 

C1FA  =  ■  ?  {^  -  1}  =  "  °*382  (7.7.1) 


c1WA  =  -rlAr^l  =  -  0-222  (7.7.2) 


'IMA  ■■  ■  o-\      ^ 


At  point  C   - 


C1FC  =  -  ?  ^  +  1)   =  -  1.U3  (7.7.3) 


ClPiC  7  +  $  \£jr\  =  +  °'825  (7.7-W 

Consequently,  the  resulting  solution  over  the  transition  range  as  represented 
by  Eq.  (7»*0  can  now  be  reduced  to  the  final  form 


^vl-5  ct'*:i) 
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which  applies  over  the  range 

-  £  (J5  +  1)  <;  C1F  ^  -  J  (vS  -  1)  (7.8.2) 

8.   Propagation  of  Vorticity  and  Development  of  the  Vortex  Ring  State 

Consider  conditions  at  station  s  far  downstream  of  the  actuator.  The  slip- 
stream at  this  location  is  moving  with  uniform  velocity  Va ,  the  ambient  flow 
with  uniform  axial  velocity  V.  In  the  idealized  inviscid  case,  the  two  irro- 
tational  flow  regions  are  separated  by  a  vortex  sheet  of  infinitesimal  thick- 
ness in  the  form  of  a  cylindrical  surface.  The  vortex  lines  are  circumferential; 
the  surface  consists  of  a  continuous  distribution  of  vortex  rings.  According 
to  the  Helmholtz  laws  of  vortex  motion  for  an  inviscid  fluid,  the  circulation 
moves  with  the  fluid  element  which  exhibits  it.  These  rotating  fluid  elements 
which  constitute  the  interfacial  vortex  sheet  themselves  move  with  a  velocity 
V  which  is  the  algebraic  mean  of  the  velocities  V  and  V_  on  either  side  of 
the  interface,  thus 

v(U  =  i(v  +  vs)  (8.1.1) 

By  comparing  this  with  Eq.  (5-8.2)  we  note  that 

V  =  V,  (8.1.2) 

uj    d 

In  other  words  the  velocity  of  propagation  V  of  the  remote  vorticity  is  the 
same  as  the  flow  velocity  V,  through  the  disc  itself!  This  applies  also  to 
the  transition  mode  T  if  we  consider  the  inner  and  outer  streams  individually. 

We  may  view  the  vortex  rings  as  vorticity  shed  originally  from  the  actuator 
disc  and  propagated  at  an  eventual  axial  velocity  V  =  V,.  Normally,  the  pro-  . 
pagation  is  in  the  sense  leading  away  from  the  actuator,  which  tends  to  produce 
a  stable  flow  pattern.   If  the  sense  of  propagation  should  become  reversed  for 

66 


any  reason  so  that  the  vortex  rings  were  to  move  toward  the  actuator  disc, 
the  tendency  would  be  to  disrupt  the  original  flow  pattern.  All  of  the 
vorticity  would  tend  to  collect  and  become  fixed  at  the  actuator,  thus 
leading  to  the  establishment  of  the  vortex  ring  regime  as  shown  earlier  in 
the  schematic  diagram  Fig.  2.6.3. 

This  idea  can  easily  be  reduced  to  explicit  mathematical  terms.  Consider 
firstly  flows  like  modes  P  and  W  in  which  the  slipstream  continues  in  the  same 
sense  as  the  undisturbed  flow.   In  this  case 


-f=  (1  +  X)  >0  (8.2.1) 


and  stability  in  the  above  sense  requires  that 


^=(l  +  |)>0  (8.2.2) 


Next  consider  flows  like  mode  R  in  which  the  slipstream  is  reversed  with  respect 
to  the  undisturbed  flow.   In  this  case 


■f--  (1  +  X)  <0  (8.2.3) 


and  stability  requires  that 


T-  f1+l)  <0  {8-2'k) 


The  above  conditions  may  be  consolidated  into  the  single  statement  that, 
for  stability,  it  is  necessary  that 
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r-rr  >0  (8-3-d 

It  can  be  seen  that  the  stability  constraint  (8.3.I)  is  satisfied  for  all 
values  of  \   except  those  lying  within  the  range 

-2  £  \  £  -1  (8.3.2) 

But  this  is  precisely  the  range  of  the  mode  U  which  was  previously  identified 
as  the  unstable  branch  of  the  formal  solution! 

Of  course,  the  two  end  points  identified  in  Eq.  (8.3.2)  are  points  of 
neutral  stability.  We  identify  these  as  points  V  and  N,  respectively.  Note 
that 

At  Point  N  \  =  -1  (8.1+.1) 

°lF=?.(1+!)*=-J  <8^2> 

ciPi  =  ?(1  +  !)2*  =  -£  .  <8-u-3> 

At  Point  V  X  =  -2  (8. 5.1) 

c1F  =  0  (8.5.2) 


c1Pi  =  0  (8.5.3) 


It  is  of  interest  to  compare  this  latter  point  with  another  significant 
operating  condition,  namely,  point  0.  Thus 
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At  Point  0  \  =  0  (8.6.1) 


C^  =  0  (8.6.2) 


Clpi  -  0  (8.6.3) 


This  brings  out  the  very  interesting  fact  that  although  points  0  and  V 
are  both  located  at  the  origin  in  the  CL.^  versus  C.   plane,  they  represent 
quite  distinct  operating  conditions.   In  particular,  point  0  is  a  state  of 
positive  stability  while  point  V  is  a  state  of  neutral  stability.   In  fact, 
point  V  represents  the  vortex  ring  state  previously  mentioned.  A  rotorcraft 
which  is  inadvertantly  brought  to  this  operating  point  suddenly  suffers  a 
catastropic  loss  of  lift! 

Consider  in  this  light  the  operation  of  an  actuator,  initially  in  mode 
R,  at  gradually  decreasing  levels  of  power  coefficient  C   . .  Assume  that  when 
the  actuator  reaches  point  C ,  operation  is  not  shifted  along  the  optimum 
transition  path  CBA,  but  is  controlled  so  as  to  continue  along  in  mode  R 
toward  point  V.   It  can  be  seen  that  as  the  operating  point  moves  closer  to 
point  V,  the  margin  of  stability  in  the  above  sense  becomes  less  until  at 
point  V  itself,  the  point  of  incipient  instability  is  finally  reached.   It  is 
evident  that  for  operation  unduly  close  to  this  point,,  a  quite  small  pertur- 
bation may  suffice  to  disrupt  the  precarious  stability,  and  lead  to  the 
vortex  ring  state. 

Moreover,  if  the  attempt  were  nevertheless  made  to  continue  the  above 
gradual  changes  further  into  the  unstable  region  VIM,  the  resulting  unstable 
flow  pattern  would  inevitably  collapse,  to  be  succeeded  in  most  cases  by  the 
vortex  ring  pattern. 
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In  contrast  with  this,  during  a  shift  along  the  optimum  transition  path 
CBA,  operation  too  near  the  points  of  incipient  instability,  V  and  N,  is 
avoided.   Hence  path  CBA  is  theoretically  stable  in  this  sense  as  well  as 
being  the  transition  path  of  highest  performance. 

Of  course,  for  a  real  actuator,  there  may  be  various  other  types  of 
instabilities  to  consider.  Fortunately,  these  are  usually  of  a  less  drastic 
kind.  Even  for  the  idealized  actuator,  however,  there  are  other  types  of 
instabilities  as  well.  One  particular  type  is  analyzed  in  section  10  of 
this  paper.   In  general,  however,  the  phenomena  of  hydrodynamic  instability 
ir  real  systems  are  extremely  complex  and  cannot  be  analyzed  by  the  elementary 
methods  of  this  paper. 

9-  Rotation  Effects  on  Actuator  Performance 
9-1  The  Idealized  Finite  Torque  Actuator 

The  ideal  actuator  considered  so  far  involves  purely  axial  flow  in  the 
remote  slipstream  and  is,  in  effect,  a  zero  torque  device.   In  principle,  this 
condition  could  be  approximated  in  various  ways .   One  method  would  be  to  use  a 
pair  of  co-axial  counter-rotating  rotors  which  exert  equal  and  opposite  torques 
upon  the  flow.  Another  method  would  be  to  use  a  rotor  and  a  co-axial  stator 
which  exert  equal  and  opposite  torques.  A  third  method  would  be  to  utilize 
a  single  rotor  at  rotational  speeds  so  high  that  the  necessary  power  could  be 
transmitted  at  torque  levels  low  enough  to  be  negligible. 

In  most  cases,  the  above  methods  are  for  one  reason  or  another  impractical, 
and  it  becomes  useful  to  consider  in  somewhat  greater  detail  the  optimum  per- 
formance which  can  be  attained  by  a  single  rotation  actuator  which  transmits 
appreciable  shaft  power  at  limited  rotational  speed  and  therefore  with  signi- 
ficant torque  effects .   In  general  such  a  device ,  which  we  term  here  the 
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finite  torque  actuator,  even  if  idealized  in  every  other  respect,  will  always 
perform  somewhatless  effectively  than  the  comparable  idealized  zero  torque 
actuator . 

9.2  Idealized  Rotation  Power  Loss 

Consider  three  actuators  all  of  which  operate  with  the  same  values  of  p, 

D,  V  and  F,  respectively,  hut  with  different  values  of  shaft  power  P.  The 

first  is  the  ideal  zero  torque  actuator  with  ideal  shaft  power  P.  as  previously 

analyzed.  The  second  is  an  idealized  finite  torque  actuator  and  the  third  is 

the  real  actuator;  both  of  these  latter  devices  operate  at  the  same  shaft  speed 

N.   The  shaft  power  P!  of  the  idealized  finite  torque  actuator  equals  (P.  +  P  ) 

where  P.,,  which  is  always  non-negative,  is  here  termed  the  idealized  rotation 

loss.  The  shaft  power  P  of  the  real  actuator  was  previously  written  as  equal 

to  (P.  +  P_)  where  P„,  which  is  always  non-negative ,  denotes  the  total  power 

loss.  This  same  total  shaft  power  P  can  now  be  represented  in  somewhat  more 

detail  by  the  expression  (P.  +  P„T  +  P  )  where  P  ,  which  we  term  here  the 

1    N    r        r 

residual  loss,  is  again  a  non-negative  quantity.  This  analysis  resolves  the 

total  power  loss  P_  into  an  idealized  rotation  loss  P„  plus  a  residual  loss 

P  . 
r 

9.3  Conventional  Coefficients  in  Terms  of  Slipstream  Parameters 

In  the  preceding  discussion  of  power  losses,  it  is  necessary  to  distinguish 
between  various  power  components  and  in  this  context  the  total  shaft  power  of 
the  idealized  finite  torque  actuator  is  denoted  by  the  somewhat  elaborate 
symbol  P.'.  However,  in  this  and  subsequent  subsections,  the  discussion  is 
restricted  solely  to  the  analysis  of  the  idealized  finite  torque  actuator  itself, 
hence  we  elect  to  denote  this  same  power  P. '  hereafter  by  the  simpler  and  more 
convenient  symbol  P.  Consequently,  the  significant  parameters  of  our  problem 
are  p,  D,  V,  F,  N  and  P.  Note  that  this  list  now  includes  N,  which  is 
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non-vanishing,  finite  and  positive.  Hence  it  becomes  advantageous  to  revert 
to  the  use  of  the  conventional  performance  coefficients  C  ,  C  and  J  based  on 
p,  N  and  D  as  dimensional  reference  parameters. 

It  is  shown  in  a  later  subsection  that  the  performance  of  the  idealized 
finite  torque  actuator  can  be  analyzed  in  terms  of  conditions  in  the  remote 
slipstream,  station  s.   It  is  necessary,  of  course,  to  distinguish  between 
the  diameter  of  the  slipstream,  which  we  here  denote  as  2R,  and  the  diameter 
D  of  the  actuator  disc  itself.   It  is  also  convenient  to  express  shaft  speed 
both  in  conventional  units  of  rev/sec  as  denoted  by  N  and  in  scientific  units 
of  rad/sec  as  denoted  by  cu.  With  this  notation,  the  required  performance 
coefficients  may  be  expressed  as  follows: 


n      -  P 


r?(^\ 


5 


'PN3D-    MD/    p*3R5 


(9.3.2) 


*■  (&  -Mt  U)  (9-3-3) 


S).k     Basic  Slipstream  Analysis 

An  actuator  which  transmits  shaft  power  P  to  an  initially  nonrotating 
stream,  and  which  itself  rotates  at  shaft  angular  velocity  u>,  necessarily 
exerts  a  torque  of  amount  P/co  upon  the  flow.  Assuming  polar  symmetry  for 
the  idealized  flow  model,  we  conclude  that  the  above  torque  sets  up  in  the 
remote  slipstream  a  distribution  of  tangential  velocity  V  which  varies  in 
some  manner  as  a  function  of  radius  r.   In  general,  the  axial  velocity  Vg 
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in  the  remote  slipstream  can  also  vary  as  some  function  of  r.  The  stream- 
surfaces  of  this  flow  are  straight  co-axial  cylinders. 

The  pressure  at  the  outer  edge  of  the  slipstream,  radius  R,  equals  the 
ambient  pressure  p.  For  a  rotating  slipstream,  the  pressure  pG  at  any  interior 
radius  r  is  below  ambient.  Considerations  of  radial  balance  of  forces  show 
that  the  amount  of  pressure  drop  is 
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'P  -  Pc 


■/* 


dr 


(9A.1) 


Applying  the  energy  equation  between  a  point  in  the  slipstream  and  a 
point  far  upstream  of  the  actuator,  we  find  the  net  energy  input  to  the  flow 
to  be 


m  = 


a  + 


v2 
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v2' 
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v2 


E  +  X- 

P    2 


(9A.2) 


The  energy  input  can  also  be  written  in  terms  of  shaft  power  input  per 
unit  mass  flow,  where  shaft  power  is  expressed  in  terms  of  torque  and  angular 
velocity,  and  where  torque  is  written  in  terms  of  change  in  angular  momentum. 
The  result  is  simply 


AE  =  (u  V  r 
u 


(9.^.3) 


Upon  combining  these  last  three  equations,  we  obtain  the  basic  relation 

R 


V 


u 


dr 


O.h.h) 


This  shows  how  the  distributions  of  V  and  V0  in  the  slipstream  must  be  related 

u     S 

for  any  assigned  values  of  the  parameters  V  and  tu. 
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9.5  Idealized  Velocity  Distribution  in  the  Slipstream 

Our  basic  aim  is  to  find  a  distribution  of  V  which  will  leave  V\, 

u  S 

unchanged  from  the  result  obtained  for  the  original  zero  torque  actuator. 
In  other  words ,  we  require  a  V  distribution  which  produces  a  uniform  axial 
velocity  V_  across  the  slipstream.  Moreover,  it  must  be  possible  to  prescribe 
the  magnitude  of  V„   arbitrarily. 

It  turns  out  that  the  required  solution  consists  of  two  distinct  regions. 
There  is  firstly  a  non-swirling  inner  core  of  radius  rn.   Outside  the  core, 
there  is  an  annular  region  of  simple  free  vortex  flow.  Thus 

For  0  £  r  £  r.    Vr  =  0  (9-5.1) 

0     u 


For  r^^r^R     Vr  =  r=  constant  (9.5-2) 

0  u 


abstituting  V„  from  Eq.  (9.5-2)  into  ($.k.k)   gives  the  result 

(9-5-3) 

With  all  other  quantities  arbitrarily  prescribed,  this  equation  fixes  the 
characteristic  swirl  constant  r  of  the  outer  free  vortex  flow. 

The  total  pressure  drop  through  the  outer  annulus  can  now  be  found  by 
substituting  the  free  vortex  relation  (9.5.2)  into  the  radial  balance  equation 
(9.U.1).   The  result  is 
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(9.5A) 


For  the  flow  through  the  nonrotating  inner  core,  there  is  no  net  energy  input 
and  the  energy  equation  gives 

7^ 


v02  -  v2  \   /p  -  p, 


(9-5.5) 


By  eliminating  the  pressure  term  between  these  last  two  expressions,  we  obtain 
a  solution  for  the  core  radius  in  the  form 


v  2  -  v2 


where  F  is  now  known  from  the  prior  solution  of  Eq.  (9.5.3). 

Eqs .  (9.5.1)  through  (9. 5. 6)  completely  define  all  features  of  the  slip- 
stream flow  which  correspond  to  arbitrarily  specified  values  of  the  parameters 
Vq ,  V  and  to. 
9.6  Slipstream  Thrust  Integral 

Applying  the  axial  momentum  law  to  the  propeller,  we  can  now  write  the 
net  thrust  in  the  form  of  an  integral  over  the  slipstream  as  follows: 


F  =  2ttp  \v  [V_  -  V  -  /  —  dr  r  dr  (9-6.1) 


The  first  term  inside  the  integrand  represents  the  basic  thrust  produced 
by  the  momentum  increase  across  the  actuator.   The  second  term,  which  can  be 
recognized  as  the  pressure  drop  integral  of  Eq.  (9.^.1) »  represents  a  loss 
of  thrust  associated  with  the  lowered  pressures  near  the  axis.   This  equation 
also  indicates  why  it  is  necessary  to  maintain  a  nonrotating  core  of  adequate 
radius  rQ.   For  if  rQ  be  allowed  to  decrease  unduly,  the  unfavorable  negative 
pressure  integral  eventually  dominates  the  solution.   The  non- swirling  core 
of  radius  rQ  eliminates  this  difficulty. 
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In  fact,  over  the  core  of  the  present  solution,  the  integrand  of  the 
thrust  equation  (9. 6.1)  reduces  as  follows 
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dr  =  vs  v  -  VI 
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>  0 


(9.6.2) 


It  is  interesting  to  observe  from  this  result  that  the  non-swirling  core  still 
makes  a  net  positive  contribution  to  the  thrust  despite  the  fact  that  it 
receives  no  energy  input  and  undergoes  no  pressure  rise  in  passing  through  the 
actuator!  The  reason,  of  course,  is  that  it  discharges  as  an  axial  jet  into 
the  region  of  reduced  pressure  created  by  the  outer  swirling  flow.   It  exerts 
this  thrust  not  as  a  pressure  i   ease  across  the  actuator,  but  as  a  suction 
force  distributed  along  the  con   -ging  stream  surface  which  separates  the 
inner  and  outer  flows .  The  outer  flow  then  transmits  this  force  to  the 
actuator  disc  itself. 
9.7  Slipstream  Power  Integral 

By  integrating  the  energy  relation  (9.*+ -3)  with  respect  to  mass  flow  rate 
across  the  entire  slipstream,  the  following  expression  is  obtained  for  shaft 
power,  namely, 
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P  =  2n  p  u>  J    Vg  V^  r2  dr 


(9.7.1) 


Likewise,  from  (9.U.I)  and  (9.U.2)  we  can  obtain  the  alternative  form 
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The  form  (9.7-1)  is  much  simpler  to  compute  ,  hut  (9.7.2)  is  interesting 

theoretically,  especially  in  conjunction  with  the  force  equation  (9. 6.1). 

Notice  that  the  effect  of  the  terms  in  V  is  to  decrease  the  thrust  delivered 

u 

and  simultaneously  to  increase  the  shaft  power  required!  The  adverse  effect 
of  slipstream  rotation  on  overall  performance  is  at  once  evident  from  these 
formulas.  Fortunately,  at  sufficiently  high  shaft  speeds,  these  effects 
become  small. 
9.8  Non-dimensional  Slipstream  Parameters 

In  applying  the  foregoing  results  it  is  useful  to  non-dimensionalize  all 
slipstream  quantities  and  relations  initially  in  terms  of  the  three  dimensional 
reference  parameters  p,  u)  and  E.   In  this  connection  we  introduce  the  following 
dimensionless  quantities,  namely, 

(9.8.1) 


(9.8.2) 
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9.9  Slipstream  Contraction  or  Expansion 

A  glance  at  Eqs .  (9*3)  now  shows  that  to  convert  the  slipstream  force  and 
power  coefficients  into  the   conventional  force  and  power  coefficients  C. 


and  Cp,  and  to  convert  slipstream  advance  ratio  h  to  the  conventional  advance 


F 
2 


ratio  J,  it  is  necessary  to  know  the  slipstream  area  contraction  ratio  (— I 
Applying  the  continuity  condition  between  actuator  disc  and  slipstream 
gives  the  general  result 

D/2 

/ 

(9-9.1) 
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However,  for  the  specific  solution  developed  in  this  section,  the  addition 

of  the  peripheral  velocity  components  V  in  the  slipstream  does  not  affect  the 

original  axial  velocity  V  .      In  fact  the  solution  was  specifically  tailored  to 

meet  this  requirement.   For  incompressible  flows  having  polar  symmetry,  a  free 

vortex  motion  superimposed  on  an  arbitrary  axial  motion  does  not  affect  the 

latter.  This  principle  is  well  known  in  fluid  mechanics  and  is  simply  utilized 

here  without  detailed  derivation  or  proof.  We  merely  note  that  the  superimposed 

free  vortex  motion  does  not  affect  the  axial  velocity  distribution  V,.  Also, 

since  Va  and  V  both  happen  to  be  constants,  Eq.  (9.9.l)simplifies  to 
o      d 

f)  'vTvT-l  (9.9.2) 

Upon  non-dimensionalizing  this  in  the  p,  tu,  r  system  we  obtain  finally 

hv\2       K  +  p 

(9.9.3) 
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9-10  Swirl  Intensity,  Loading  Limit  and  Net  Operating  Range 

Further  details  of  the  non-dimensionalized  solution  now  work  out  as 
follows.  The  basic  relation  (9-5.3)  which  fixes  slipstream  swirl  intensity 
takes  the  form 


,  +  £)  A  =  Y  -  i  Y* 


(9.10.1) 


It  is  useful  to  solve  this  for  y  "by  the  quadratic  formula,  and  to  arrange 
the  resulting  solution  in  the  following  special  form,  namely, 


i  +  |  A 


(9-10.2) 


where  the  auxiliary  parameter  Q   is  defined  by  the  expression 


(9-10.3) 


It  is  also  convenient  for  some  purposes  to  expand  this  solution  into 
series  form  as  follows. 


c  =  i  +  i  [h  +  |)  a  +  i  ft  + 1)2  a2  + 


(9-io.U) 


Inspection  of  Eq.  (9. 10. 3)  discloses  an  important  physical  operating  limit, 
Since  £  is  always  real,  the  quantity  under  the  radical  can  never  become  negative, 
We  conclude  that 
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(9.10.5) 
(9.10.6) 

(9.10.7) 
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Solving  Eq.  (9. 10. 5)  for  the  loading  limit  iV^  by  the  quadratic  formula 
gives  the  result 


=  +  J 1  +   k  -  K 


Amax=+V1+h  -  h  (9.10.8) 

There  is  also  a  negative  root  for  A  which  would  be  significant  if  we 
were  considering  reversed  flows,  mode  R,  in  this  analysis.  However,  the 
present  discussion  is  limited  to  the  simple  modes  P  and  W  only,  and  in  this 
context  the  negative  root  is  not  relevant. 

On  the  other  hand  there  is  another  restriction  on  A  arising  from  the  very 
fact  that  the  analysis  of  this  section  is_  specifically  limited  to  modes  P  and 
W.   In  terms  of  the  earlier  parameter  X,  we  may  say  that  the  present  study  is 
restricted  to  the  range 

.2  f^i  :  M  *  x  £  +  *  (9.10.9) 

V  S 


where 


x  -  -f  =  \  (9.10.10) 


Upon  combining  these  conditions,  we  find  that  the  present  analysis  is 
restricted  to  the  net  operating  range 


_  2     V- 2  -  1)  H*AsC  +  </l  +  K-H  (9.10.11) 


where  positive  values  of  A  correspond  to  mode  P,  negative  values  to  mode  W. 
9.11  Slipstream  Core  Radius 

For  any  arbitrarily  assigned  values  of  h  and  A  which  satisfy  the  constraint 
(9.10.11),  parameter  £  is  fixed  by  (9. 10. 3)  or  (9.10.U)  and  parameter  y   by 
(9.10.2).  Eq.  (9.5.6)  which  fixes  the  relative  radius  T^  of  the  non-swirling 
core  in  the  slipstream  can  then  be  non-dimensionalized  in  the  form 
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(9.11.D 


9.12  Parametric  Solution  for  Conventional  Performance  Coefficients 

With  all  of  the  foregoing  details  fixed,  it  is  a  routine  matter  to  non- 
dimensionalize  and  integrate  the  fundamental  thrust  and  power  relations  (9-6.1) 
and  (9.7.1).  Upon  utilizing  relations  (9«3)»  we  finally  obtain  the  following 
parametric  solutions  for  C„,  C  and  J. 
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(9-12.1) 


(h  +  A)  Y  1  -  V 


(9-12.2) 


J  =  TT  H 


A 
H  +  A 


(9.12.3) 


As  explained  above,  the  two  auxiliary  parameters  y  an^  \>  which  occur 
in  these  equations  are  themselves  known  functions  of  the  two  independent 
parameters  h  and  A.  Hence  in  principle  these  equations  suffice  to  fix  a  family 

of  curves  of  C„   versus  J,  with  C^  as  the  parameter  which  changes  from  one  curve 

F  P 

to  the  next. 

9.13  Static  Thrust  and  Power  Coefficients 

The  reader  can  easily  verify  that  under  static  conditions  (V  =  0),  the 
foregoing  parametric  solution  reduces  to  the  following  form. 
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(9.13.6) 


Values  of  C„  ,  Cp  and  C^/Cp    as  calculated  from  these  equations  were 
previously  presented  in  Fig.  2.9.I 
9-lU  Maximum  Attainable  Thrust  and  Power  Coefficients 

The  reader  can  also  readily  verify  that  under  the  absolute  maximum  loading 
limit  of  Eq.  (9.10-8)  the  general  parametric  solution  reduces  to  the  following 
simplified  results. 
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Eq.  (9« l^-3)  shows  that  under  these  absolute  maximum  loading  conditions 

the  non-swirling  core  occupies  exactly  half  the  cross -sectional  area  of  the 

slipstream!  Under  normal  light  loadings,  the  core  is  of  course  very  much 

smaller  than  this . 

Curves  of  C.J    and  C  ]    versus  J  as  computed  from  these  relations 
/MAX       /MAX 
were  previously  presented  in  Fig.  2.9.2. 

9.15  Relative  Flow  Angle  Downstream  of  Actuator 

In  studies  relating  to  pitch  setting  of  variable  pitch  actuators ,  it  is 

frequently  useful  to  be  able  to  compute  the  theoretical  relative  flow  angle 

B,  at  radius  r,  on  the  downstream  side  of  the  actuator  disc.  For  the  idealized 
d  d 

actuator,  this  is  defined  by  the  expression 


V. 


tan  B,  = 


d   r 


(9-15.1) 


To  relate  this  to  corresponding  values  in  the  remote  slipstream,  we  note 


that 
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rd  = 
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Upon  combining,   rearranging  and  non-dimensionalizing  we  finally  obtain 


tan   3. 
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The  corresponding  flow  angle  (3  in  the  remote  slipstream  itself  is  given 
by  the  simpler  expression 


tan  (3  = 


1  -X 
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(9.15.5) 


Both  of  these  angles  vary  as  functions  of  Tj  .  To  find  an  angle  which  is 
respresentative  of  the  pitch  setting  in  some  approximate  overall  sense,  it  is 
necessary  to  use  an  appropriate  mean  relative  radius  T)  .  For  example,  we  may 
choose  as  a  suitable  mean  value  for  this  purpose  the  quantity 


rl  + 


^=T\n 
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(9.15.6) 


For  an  actuator  with  a  sufficiently  large  number  of  blades  (high  blade 
density)  the  outlet  flow  angle  R,  at  radius  TL  will  approximately  equal  the 
outlet  blade  angle  itself  at  this  point.  For  larger  blade  spacings,  the 
relation  between  flow  angle  and  blade  angle  is  more  complex  and  lies  outside 
the  scope  of  this  discussion. 
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By  combining  relations  (9.15.6)  with  (9.15.U)  and  imposing  the  constraint 
B_  =  constant  (at  radius  71  )  (9.15.7) 

it  is  possible  to  find  a  relation  between  h  and  A  which  corresponds  approximately 
to  constant  pitch  setting  (at  high  blade  density).  The  corresponding  values 
substituted  into  the  general  parametric  solution  of  section  9«12  then  roughly 
simulates  the  performance  of  an  actuator  at  fixed  pitch  setting. 
10.  Dynamic  Stability  of  the  Idealized  Lifting  Rotor  in  Vertical  Ascent  or 

Descent 

In  this  section  we  consider  the  response  of  a  lifting  rotor  to  small 
perturbations  from  equilibrium.  The  treatment  is  deliberately  simplified  to 
forcus  attention  strictly  on  the  basic  phenomena.  An  ideal  actuator  is  assumed, 
with  P„  =  0.  We  analyze  a  rotorcraft  of  gross  weight  W  operating  in  steady 
vertical  climb  or  descent,  except  for  the  imposed  perturbations.  At  the  steady 
reference  condition,  lift  L  is  equal  to  weight  W.  But  unlike  the  weight  ¥ 
which  is  constant,  lift  L  is  subject  to  small  perturbations  6L  about  the 
equilibrium  value.   Lift  L  and  weight  ¥  are  the  only  forces  considered  in  this 
analysis,  drag  and  other  incidental  forces  being  deliberately  ignored. 

Under  these  circumstances  it  is  advantageous  to  substitute  ¥  for  L  as 
dimensional  reference  parameter.  Hence  we  now  work  in  the  p,  D,  ¥  dimensional 
system.  The  relevant  coefficients  are  defined  as  follows 

C   =  — 2—  (10.1.1) 

1 

D 


C™  "  -V  (10.1.2) 

3H   1  /w 


D  V  p 
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C3L=f  (10'1-3) 


In  this  system,  the  fundamental  equations  for  the  ideal  actuator  take 
the  following  form  - 


For  Modes  P,  R  and  W  - 


C„  C  •       7CL,  C. 


%-^*WM  +  ^  (M'-2;l) 


For  Mode  T 


°3P  =  •  335  T  ^  "  °3L  C3H  (10-2-2) 

For  the  stability  analysis  it  is  necessary  to  find  the  partial  derivatives 

ac   \        /ac  \ 

—  and  y,0  .   This  is  easily  accomplished  by  differentiating  Eqs.  (10.2) 


\  ^36/     \  ^31 
However,  in  the  end  we  require  only  the  equilibrium  values  of  these  various 

quantities.  The  equilibrium  values  are  obtained  simply  by  setting  C_  =  1 

in  the  various  expressions.  Of  course,  this  may  be  done  only  after  the 

required  differentiations  are  carried  out.   In  this  way  the  following  results 

are  obtained 

For  modes  P,  W  and  R  - 


(10.3.1) 
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ac 


y  = 


3P   H 


_H 

LV  2  > 


sc 


TT  Cp.. 


3L 


4  «.** 

2  /     it  cp 


(10.3.3) 


For  mode  T  - 


C3P  =  "  4  T  CH3  •  °H 


(lO.U.i) 


X  =^3P=  -  -HTc.2  -  1 
^H    n#   H 


(10. U. 2) 


Y  "  *3L     H 


(10. U. 3) 


The  basic  equations  (10.2)  can  be  represented  symbolically  in  the  form 


V  C3l)  (10-5-1) 


Consequently,  the  small  changes  occurring  during  a  perturbation  may  be  written 


ac 


6C 


3P 


8C. 


W^3P 


6C. 


\  dC-A  /   3H   \BC«,y   3L 


or 


6C3p  =  X  6C3-  +  Y  5C3L 


(10.5.2) 


The  equation  of  motion  of  the  rotorcraft  during  a  small  perturbation  is  - 


_6L  =  1  _d_ 
W    g  dt 


(8H) 


(10.6.1) 
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By  utilizing  the  definitions  (10.1.2)  and  (10. 1. 3)  this  may  be  put  into  the 
form 


6C 


_  1  _d_ 
3L  "  g  dt 


D  .  p  ^H 


(10.6.2) 


It  is  now  advantageous  to  introduce  a  non-dimensional  time  9  defined  as  follows 


0  = 


_1_  /W 
gD  J  P 


(10.6.3) 


When  expressed  in  terms  of  this  dimensionless  time  coordinate  9,  the  equation 
of  motion  (10.6.2)  becomes  simply 


6C3L  "  d6 


6C3H 


(10.6. h) 


It  is  instructive  to  include  in  the  analysis  the  possibility  of  incorporating 
an  automatic  power  control  system.  Specifically,  we  assume  that  power  input 
or  output  may  be  modulated  in  response  to  perturbation  velocity  and/or  accel- 
eration as  sensed  by  appropriate  instrumentation.  Thus  let 


6C3P  =  ■  *  d9 


6C 


3H 


-  v 


6C 


3H 


(10.7.1) 


where  p,  and  v  are  two  non-negative  constants  representing  the  response  of  the 
control  system  to  perturbation  acceleration  and  perturbation  velocity, 
respectively.  Of  course,  we  may  set  jj,  =  v  =  0,  thus  simulating  a  control 
system  which  maintains  constant  power  independently  of  the  state  of  the  motion 
of  the  craft. 

Upon  substituting  Eqs.  (10. 6. k)   and  (10. 7.1)  into  (10.5.2),  we  obtain  the 
basic  response  equation 
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de 


6C 


3H 


X  +  y 

Y  +  m, 


6C  . 
3H 


(10.8.1) 


The  solution  of  Eq.  (9-8.1)  can  "be  taken  in  the  form 


6C 


3H 


6C 


3H 


36 


0 


(10.8.2) 


_d_ 

de 


6C 


3H 


=  0 


SC 


3H 


.pe 


^0 


(10.8.3) 


Upon  substituting  these  last  two  expressions  into  the  basic  response 
equation  (10.8.2),  we  finally  obtain  the  result 


3  =  - 


X  +  v, 
Y  +  u. 


(10. 8. k) 


where  the  functions  X  and  Y  are  as  defined  in  Eqs.  (10. 3)  and  (10. U),  while  \i 
and  v  are  two  non-negative  constants  which  characterize  the  control  system, 
if  any.  Under  conditions  of  constant  power  input  or  output,  p,  =  v  =  0,  and 
Eq.  (10.8.U)  reduces  further  to 


3  =  - 


x 


(10.8.5) 


Stability  of  operation  in  vertical  ascent  or  descent  requires  that  3 
be  negative.  We  shall  say  that  Eq.  (10.8.5)  defines  the  inherent  stability 
of  the  idealized  rotor,  while  (10. 8. k)   defines  the  augmented  stability. 

Study  of  the  actuator  derivative  functions  X  and  Y  as  given  earlier  in 

Fig.  2.10.1  shows  that  X  is  positive  for  all  positive  and  negative  values  of 

the  vertical  velocity  Corr.  The  function  Y,  on  the  other  hand,  turns  out  to 

3H 


89 


be  positive  for  values  of  C  •  algebraically  greater  than  — —  =  -  1.693,  and 

3H  -   JZ  ' 

negative  for  values  of  C0«  algebraically  less  than  — —  .  This  cross -over 
point  coincides  with  point  M  on  the  operating  envelope. 

We  conclude  that  the  ideal  lifting  rotor  is  inherently  stable  over  the 
range 

-  4s  ^  C  •  ^  +  oo  (10.8.6) 

This  includes  all  rates  of  ascent,  as  well  as  descent  over  the  full  range  of 
power-on  reversed  thrust  operation,  and  descent  over  an  initial  small  portion 
of  the  windrailling  range.   It  is  to  be  expected  that  for  any  real  rotor,  the 
range  of  inherent  stability  is  likely  to  be  considerably  more  restricted  than 
this. 

Eq.  (10. 8. k)   shows  very  clearly  that  the  range  of  stable  operation  can  be 
extended  as  far  below  the  operating  point  M  as  desired  simply  by  introducing 
the  requisite  amount  of  acceleration  response  as  denoted  by  the  constant  (i. 
Note  that  the  velocity  response  factor  v  does  not  affect  the  range  of  stability 
but  merely  increases  the  rate  at  which  disturbances  either  grow  or  decay. 

Typical  stability  response  curves  for  various  combinations  of  control 
parameters  as  computed  from  the  foregoing  equations  were  shown  earlier  in 
Fig.  2.10.2. 
11.  Actuator  Efficiency  and  Effectiveness 

In  analyzing  the  performance  of  a  propeller,  windmill  or  lifting  rotor, 
we  are  concerned  not  only  with  determining  the  absolute  values  of  force,  power 
and  so  on,  but  also  with  assessing  the  comparative  performance  of  the  device 
in  relation  to  its  theoretical  potential.   In  this  connection,  an  essential 
first  step  is  to  reduce  all  performance  parameters  to  appropriate  dimensionless 
form  since  this  normalizes  the  performance  of  devices  of  every  size  and  capacity 
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to  a  common  basis  of  comparison.  However,  this  first  step,  while  essential, 
is  not  necessarily  sufficient  in  itself.  For  example,  to  know  the  actual 
values  of  the  conventional  coefficients  C  ,  C  and  J  in  a  specific  instance 
is  indispensable  for  some  purposes,  but  nevertheless  does  not  directly  indicate 
the  degree  of  departure  from  ideal  performance.  Fortunately,  this  significant 
aspect  of  actuator  performance  can  be  expressed  simply  and  consistently  by  the 
concept  of  power  loss  P-  as  defined  elsewhere  in  this  paper.  Naturally,  this 
loss  itself  should  be  expressed  in  dimensionless  form  by  the  coefficient  C„, 
Clf  C2f  or  C3f  as  ^Propri3^6 »  according  to  the  context. 

In  some  cases  it  is  useful  to  break  the  total  shaft  power  P  down  in  some- 
what more  detail  as  shown,  for  example,  in  Table  11.1.  From  this  breakdown  of 
total  shaft  power  P  into  the  four  basic  components  shown,  it  is  possible  to 
define  a  sequence  of  significant  power  ratios  as  indicated  by  Eqs.  (ll.l.l) 
through  (11.2.3)  of  the  table.  Note  that  these  definitions  all  employ  the 
sign  changing  operator  <j.   whose  definition  is  repeated  in  Eq.  (ll.U.l).  The 
purpose,  of  course,  is  to  keep  the  efficiency  and  effectiveness  ratios  as 
defined  in  the  table  always  less  than  unity,  regardless  of  the  algebraic  sign 
of  the  ideal  shaft  power  P.  . 

The  key  relations  summarized  in  Eqs.  (11. 3)  follow  at  once  from  the  earlier 
definitions.  Notice  from  Eq.  (11. 3. 3)  that  the  overall  power  effectiveness  ep 
is  closely  related  to  the  power  loss  coefficient  Cp-p-. 

The  familiar  concept  of  efficiency,  when  applied  to  actuators,  is  very 
useful  in  some  contexts,  but  completely  useless  in  others.  For  example,  for 
a  propeller  operating  under  cruise  conditions,  overall  efficiency  Tl  is  ordinarily 
a  factor  of  prime  importance.  On  the  other  hand,  for  a  lifting  rotor  in 
stationary  hovering  flight,  both  T|  and  TL  vanish  identically  and  obviously  fail 
to  provide  relevant  figures  of  merit.  Nevertheless,  the  power  effectiveness 
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factors  e  ,  e  and  e  continue  to  be  significant  in  this  case  and  do  provide 
useful  figures  of  merit. 

Notice  that  the  definition  of  overall  power  effectiveness  ep  amounts  to 
a  comparison  of  ideal  and  actual  power  under  conditions  of  equal  thrust.  In| 
fact  the  definition  (11.2.3)  may  be  rewritten  in  equivalent  form  as 


eP  = 


-.  ai 


'lPi 


'IP 


(11.5) 


J  C 


IF 


where  the  subscript  C.^  signifies  that  ideal  and  actual  power  coefficients  ai 

XT 

compared  on  the  basis  of  equal  force  coefficient. 

It  is  equally  possible,  however,  to  compare  the  ideal  and  actual  force 
coefficients  on  the  basis  of  equal  power  coefficients  for  both,  and  in  some 
situations  this  is  a  preferable  method  of  comparison.  We  recall  that  for 
any  assigned  value  of  C.  p,  the  ideal  actuator  has  two  distinct  values  of 


thrust  coefficients  which  we  here  denote  as  C 


IF 


and  C 


MAX 


IF) 

/MTTJ 


At  the  same 


value  of  C_ _,  the  actual  thrust  coefficient  C.^   always  his  within  the  range 


w^^H 


(11.6.3 


Consequently,  there  are  two  significant  force  ratios  which  we  define  as  folic 
namely , 


and 


C1F  1 

CF  " 

cif) 

/MAX 

"  C1F  1 

CB  " 

cif| 
Imin 

(11.6. 2 


Jc 


IP 


(11.6.3 


'IP 
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Notice  that  the  sign  changing  exponent  a  is  needed  in  Eq.  (11.6.2)  to  keep  e_ 
always  less  than  unity  regardless  of  the  sign  of  C.  p,  but  that  no  sign  changing 
exponent  is  needed  in  Eq.  (11. 6. 3).   In  this  connection,  recall  that  C.^N 

has  the  same  sign  as  C__,,  but  that  C_J\    is  always  negative.  We  may  term 

1P  "/MOT 

e  the  thrust  effectiveness  for  positive  values  of  C   ,  or  the  drag  effective- 
ness for  negative  values  of  C__.  We  may  term  e_,  the  braking  effectiveness  in 
all  cases.  This  terminology,  while  arbitrary,  is  reasonably  convenient  and 

clear.  The  values  of  e„  or  e  can  be  readily  inferred  in  any  given  case  from 

r      .D 

the  basic  relations  presented  elsewhere  in  this  paper. 

The  power  loss  coefficients,  the  efficiencies,  or  the  power  and  force 
effectiveness  ratios  defined  in  this  section  provide  suitable  figures  of  merit 
for  virtually  every  situation. 
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Table  11.1  Significant  Power  Components  and  Ratios 


Components  of  Total  Shaft  Power 

Component 

Sign 

Symbol 

Subtotals 

Thrust  Power 

+ 

FV 

1   '■ 

Basic  Slipstream  Loss 

+ 

ps  ■ 

) 

Idealized  Rotation  Loss 

+ 

'"  } 

1       * 

Residual  Loss 

+ 

',  J 

Total  Shaft  Power 


Significant  Power  Ratios 
Tl  = 


a. 

r  -.  1 


FV 
P 


\- 


eN  = 


FV 

P. 

l 


P  +  P 
i    N 


— 

-i 

e  = 
r 

l     N 

P 

eP  = 

"P." 

1 

P 

0". 

1 

Key  Relations 


Overall  Efficiency 


Ideal  Efficiency 


Tl  =  \ 


eP  =  eN  er 


[■ 


-.a. 


ep=     1  -  oC2pf 


(Table  continues  on  next  page) 
9k 


(ll.l.l) 


(11.1.2) 


Rotation  Power  Effectiveness       (11.2.1) 


Residual  Power  Effectiveness       (11.2.2) 


Overall  Power  Effectiveness        (11.2.3) 


(11.3.1) 
(11.3.2) 

(11.3-3) 


Table  11.1  Cont'd 


Sign  Changing  Operators 


P. 
a.=-jp^=  +  l  (ll.U.l) 


a  =  -fit  =  +  1  (11.  k.  2) 


w - 


95 
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This  paper  deals  with  certain  classical  fundamentals  of  fluid  mechanics 
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as  expounded  in  the  following  paper.  A  further  bibliography  on  dimensional 
analysis  may  also  be  found  therein. 
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